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ABSTRACT 

We present an analysis of the ultraluminous X-ray source (ULX) population in 75 Virgo cluster 
late-type galaxies, including all those with a star formation rate z І Мо уг”! and a representative 
sample of the less star-forming ones. This study is based on 110 observations obtained over 20 
years with the Chandra X-ray Observatory Advanced Camera for Imaging Spectroscopy. As 
part of a Large Chandra Program, new observations were obtained for 52 of these 75 galaxies. 
The data are complete to a sensitivity of 10% erg 571, with a typical detection limit of ~3 
х10%% erg 87! for the majority of the sources. The catalogue contains about 80 ULXs (0.3— 
10 keV luminosity >10° erg 571), and provides their location, observed flux, de-absorbed 
luminosity, and (for the 25 most luminous ones) simple X-ray spectral properties. We discuss 
the ULX luminosity function in relation to the mass and star formation rate of the sample 
galaxies. We show that recent models of low-mass plus high-mass X-ray binary populations 
(scaling with stellar mass and star formation rate, respectively) are mostly consistent with our 
observational results. We tentatively identify the most luminous X-ray source in the sample (a 
source іп IC 3322A with Lx = 6 x 10 erg 87!) as a recent supernova or its young remnant. 
The properties of the sample galaxies (morphologies, stellar masses, star formation rates, total 
X-ray luminosities from their point-source population) are also summarised. 


Key words: galaxies: clusters: individual: Virgo — galaxies: spirals — X-rays: binaries — X-rays: 
galaxies 


1 INTRODUCTION Lehmer et al. 2019). On the one hand, accreting compact objects 
probe the last stages of stellar evolution and collapse. On the other 
hand, binary systems with compact objects include the progenitors 
of gravitational wave mergers. Thus, studying the number and distri- 
bution of X-ray binaries provides empirical constraints to population 
synthesis models and theoretical predictions of binary merger rates 
(Belczynski et al. 2016; Kruckow et al. 2018; Spera et al. 2019). 


The population properties of compact X-ray sources are a tracer 
of star formation rate (SFR) and stellar mass (M+) (Grimm et al. 
2003; Gilfanov 2004; Mineo et al. 2012a; Fragos et al. 2013a; 
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Reliable scaling relations between galaxy properties and the total lu- 
minosity of the X-ray binary population are also essential to estimate 
the relative contribution of off-nuclear compact objects (compared 
with massive stars and quasars) to cosmic reionization and to the 
heating of the intergalactic and interstellar medium at high redshift 
(Mirabel et al. 2011; Fragos et al. 2013b; Kaaret 2014; Artale et al. 
2015; Lehmer et al. 2016). 

More specifically, here we are interested in the high-luminosity 
end of the X-ray binary population, near and above the Edding- 
ton limit of typical stellar-mass accretors (Lx > 1039 erg 571). 
This class of ultraluminous X-ray sources (ULXs: see Kaaret et al. 
2017 for a review) is populated both by confirmed neutron stars 
and candidate stellar-mass black holes (BHs), but the relative con- 
tribution of the two types of accreting compact objects remains 
uncertain and widely debated (Bachetti et al. 2014; Walton et al. 
2018; Wiktorowicz et al. 2019; Mushtukov et al. 2019). ULXs are 
found both in young and in old stellar environments, although the 
younger systems are more common, especially at highest luminosi- 
ties (2 5x 1039 erg s7!) (Swartz et al. 2004; Maccarone et al. 2007; 
Zhang et al. 2012; Lehmer et al. 2019). In young-population ULXs, 
the donor stars are usually consistent with massive blue or red super- 
giants (Lau et al. 2019; Lopez et al. 2020), but other types of donors 
with strong winds are also possible: for example, a Wolf-Rayet star, 
as seen in a ULX inside the Circinus galaxy (Qiu et al. 2019), or 
a main-sequence B-type star, as in M82 X-2 (Fragos et al. 2015). 
Mass transfer from massive donors may occur via classical wind, fo- 
cused wind (“wind Roche-lobe-overflow’’), or Roche lobe overflow 
(if the mass ratio between donor and accretor is below the thresh- 
old for common-envelope formation); the relative contribution of 
the three processes is also still actively debated (Pavlovskii et al. 
2017; Artale et al. 2019; El Mellah et al. 2019; Quast et al. 2019; 
Wiktorowicz et al. 2021). Instead, in old-population ULXs, the 
donor star may be a Roche-lobe-filling low-mass or intermediate- 
mass star; such systems would be analogous to transient low-mass 
X-ray binaries. If the Roche lobe of the compact object, and therefore 
its accretion disk, is very large, the thermal-viscous disk instability 
is predicted to cause long outburst phases and super-Eddington peak 
luminosities (Lasota et al. 2015; Hameury & Lasota 2020). Alter- 
natively, even intermediate- and low-mass stars can go through evo- 
lutionary phases when they can transfer > 1076 Mo yr! onto their 
compact companion for a sustained period of time of ~10°-10° 
yr (Wiktorowicz et al. 2017; Misra et al. 2020), for example when 
they cross the Hertzsprung gap or as they ascend the red giant or 
asymptotic giant branch. Moreover, in at least one case, for an old- 
population ULX in the Virgo giant elliptical M49, the donor star is 
a white dwarf (Maccarone et al. 2007, 2010; Steele et al. 2014). 

Most of the empirical constraints to theoretical models of ULX 
formation come from X-ray surveys of accreting compact objects 
in nearby galaxies, followed by optical/IR studies of their stellar 
counterparts and environments. The Chandra X-ray Observatory 
is the best instrument available for such surveys, thanks to its still 
unsurpassed sub-arcsec spatial resolution and very low background 
noise. Already a few years after its launch, Chandra-based surveys 
showed the difference between the ULX luminosity function (LF) 
in spiral and elliptical galaxies. The early-type-galaxy LF is an in- 
dicator of М, (Kim & Fabbiano 2004; Gilfanov 2004; Zhang et al. 
2012), while the late-type-galaxy LF is mainly proportional to SER 
(Grimm et al. 2003; Mineo et al. 20122); the young-population LF 
is a power-law with an exponential cut-off at Ly ~ 2 x 10% erg s^! 
(Swartz et al. 2011; Mineo et al. 2012a). More generally, the LF of 
any given galaxy is a mix of the two components (Lehmer et al. 
2019) and depends on the specific SFR (sSFR = SFR/M.). There 


is an additional smaller dependence of the high mass X-ray bi- 
nary (HMXB) population on the metal abundance of the host 
galaxy, in the sense that metal-poor galaxies contain more luminous 
HMXBs (Prestwich et al. 2013; Brorby et al. 2014; Kovlakas et al. 
2020; Lehmer et al. 2021). 


There are now a few hundred nearby galaxies observed by 
Chandra deep enough to reach a completeness limit of about 1039 
erg s^! or below, which can be used for ULX population studies 
(for example, a sample of about 300 galaxies within 40 Mpc was 
used by Kovlakas et al. 2020 for their study of the ULX population). 
Surveys of nuclear X-ray sources suspected of being intermediate- 
mass BHs have also been conducted with archival Chandra data 
(e.g., Chilingarian et al. 2018). A caveat is that archival Chandra 
targets are biased in favour of "interesting" galaxies, sometimes with 
higher-than-average X-ray luminosities or already known to contain 
particularly luminous ULXs discovered by earlier X-ray telescopes. 
This may affect the inferred scalings of the LF with sSFR. For this 
work, we chose instead a physically selected sample of galaxies, 
representative of the spiral population in the Virgo Cluster. The 
methodology of our sample selection is described in Section 2.1. 


The Virgo Cluster is an interesting environment for a point- 
source X-ray population study in galaxies with a range of Hub- 
ble types and star-forming activity. It is the richest and biggest 
galaxy cluster in the local universe, with a total mass (gas, stars and 
dark matter) Mtot ~ 6-8 x10!4 Мо (Kashibadze et al. 2020, 2018; 
Shaya et al. 2017; Planck Collaboration 2016). It is close enough 
(d = 16.1 + 1.0 Mpc for its central galaxy, M 87, in the weighted- 
average estimate of de Grijs & Bono 2020) that Chandra imaging 
observations can reach Lx ~ a few 1038 erg s7! in <10 ks and are 
not affected by source confusion. Virgo is also just close enough to 
permit at least a detection and possibly a rough colour identification 
of the optical counterparts; at a distance modulus of 31 mag, the 
typical optical brightness of young-population ULXs (My ~ -7 
mag to —4 mag: Ambrosi & Zampieri 2018; Gladstone et al. 2013; 
Tao et al. 2011) corresponds to V - 24-27 mag. Moreover, Virgo 
already has a huge wealth of multiband datasets, which enable us 
to determine the SFR and М, of many of its galaxies. 


A sample of 100 early-type Virgo galaxies (including E, dE 
and SO morphological types) was observed by Chandra as part of 
the “Active galactic nucleus MUltiwavelength Survey of Early-Type 
Galaxies” (AMUSE) Virgo study (Gallo et al. 2008, 2010). The 
sample was selected to coincide with the Hubble Space Telescope 
(HST) Virgo Cluster Survey (Côté et al. 2004). In turns, the HST 
sample was selected from early-type galaxies with blue brightness 
Вт < 16 mag, a published central velocity dispersion, and по pecu- 
liar properties such as mergers or close interactions. The total stellar 
mass (М.) of the AMUSE-Virgo sample was M. = 6 x 102 Мо 
(derived by Gallo et al. 2008, 2010, from the go and zo luminosi- 
ties, according to the relation of Bell et al. 2003). The Chandra 
coverage was designed to be complete down to ~ 4 х 10°8 erg s~!. 
The AMUSE X-ray survey provided an exceptional view of nuclear 
BH activity and Eddington ratio in “normal” early-type galaxies 
(Gallo et al. 2010; Miller et al. 2012; Graham & Soria 2019), as 
well as a sample of bright off-nuclear sources. 


From the AMUSE data, Plotkin et al. (2014) found 55 sources 
with an apparent 0.3-10 keV luminosity Lx > 10°? erg 571, in 
the Virgo ellipticals. However, only 13 of them had an apparent 
(projected) luminosity higher than ~ 2 x 1039 erg 571, and those 
brighter sources were statistically consistent with background AGN 
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Figure 1. Sky map of the Virgo galaxies. Stars mark the positions of the 75 
spirals in our study. Black dots mark the position of the 100 AMUSE-Virgo 
early-type galaxies; the two magenta dots represent M 87 (further to the 
north) and M 49 (further to the south). Small-sized stars represent spirals 
in our sample with low stellar masses (log(M../Mo) < 9.8); intermediate- 
sized stars represent intermediate masses (9.8 < log(M./Mo) < 10.4); 
and large-sized stars represent the most massive spirals (log(M,./Mo) > 
10.4). Red-coloured stars denote galaxies with low SFR (SFR < 0.5Mo 
yr-!, according to the WISE 12-ит proxy); green-coloured stars are for 
intermediate rates (0.5Mo yr! < SFR < 1.5Мо yr ly blue-coloured 
stars are for high rates (SFR > 1.5Мо уг”). (The choice of the three bands 
in stellar mass and SFR is arbitrary, for illustration purposes.) A dotted line 
marks the boundaries of the Virgo Cluster Survey (Binggeli et al. 1984). 
The only spiral in our sample that is outside the VCC region is included in 
the Extended VCC (Kim et al. 2014) and is a probable member based on its 
redshift-independent distance and recession speed. 


projected inside the D55 area! of the Virgo galaxies. Plotkin et al. 
(2014) concluded that in the old stellar population, there is ~ 1 
ОХ per M. = 1.6x 10!! Mo. This is consistent with the expected 
LF of low mass X-ray binaries (LMXBs) in early-type galaxies 
(Zhang et al. 2012). 

To have a more complete understanding of the nuclear and 
off-nuclear X-ray sources in Virgo, we must integrate the exist- 
ing study of its elliptical galaxies with a new study of its spi- 
ral galaxies. Virgo is a spiral-rich cluster, a signature of an early 
evolutionary stage (clusters such as Coma, by contrast, are more 
evolved and spiral-poor). However, many of the Virgo Cluster’s spi- 
rals have lost HI gas via ram pressure stripping and tidal stripping, 
and as a result, the sSFR of Virgo spirals tends to be lower than 
for field spirals of similar Jeans-Lundmark-Hubble type and mass 
(see Boselli & Gavazzi 2006 for a detailed review and discussion 


! The Dos region of a galaxy is defined as the 25th magnitude B-band 
isophote. 
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Figure 2. Distribution of radial velocities with respect to the cosmic mi- 
crowave background, for the galaxies in our sample. The mean and median 
velocities, and the Gaussian standard deviation, are also labelled in the plot. 


of such issues, and Graham 2019 for a comprehensive discus- 
sion of morphological types). For the off-nuclear X-ray sources, 
Virgo spirals contain two populations of stellar-mass objects. First, 
a sample of LMXBs in spirals; the number and luminosity of such 
sources as a function of stellar mass may not be the same as for 
the LMXBs in ellipticals. Second, a sample of young high mass 
X-ray binaries (HMXBs), in the galaxies with a higher SFR. Third, 
for centrally-located, i.e., nuclear, BHs, a Chandra study can re- 
veal, for example, the difference in the Eddington ratio distribution 
between galaxies dominated by hot gas (early-type) or cold gas 
(late-type). It can also reveal the presence of candidate X-ray-bright 
intermediate-mass BHs in the nuclei of spiral galaxies with a small 
bulge or no bulge (as predicted by galaxy/BH scaling relations), 
which can then be the subject of subsequent multiband investiga- 
tions (Chilingarian et al. 2018; Mezcua et al. 2018). Furthermore, 
the unsurpassed sub-arcsec astrometric accuracy and sharp point- 
spread function of Chandra enables a search for slightly off-centre, 
“wandering” nuclear BHs (Khan et al. 2020; Ricarte et al. 2021). In 
this paper, we will mostly focus on off-nuclear sources, and we will 
leave the analysis of the nuclear BHs and candidate intermediate- 
mass BHs to separate work. 

By contrast with the early-type galaxies, the archival Chandra 
coverage of late-type Virgo galaxies was, until recently, surprisingly 
scant. Less than 30 Virgo spirals had been observed before 2017, 
an insufficient number for population studies. Thus, we planned 
a larger survey and successfully applied for 550 ks of Chandra 
time (Cycle 18 Large Program, PI R. Soria). In the next Section, 
we describe the population properties of the Chandra sample. In 
Section 3, we outline our X-ray data analysis. In Section 4, we 
present the results for the off-nuclear source population, and place 
it in the context of the properties of their host galaxies. A census 
of the nuclear sources, a study of the diffuse X-ray emission in the 
most star-forming galaxies, and an investigation of the multiband 
counterparts will be presented in further work. 


2 A CHANDRA SAMPLE OF VIRGO SPIRALS 
2.1 Selection criteria 


The selection of targets for our study resulted from a necessary 
compromise between various ideal requirements: a large number 
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of galaxies, particularly with high SFR (because that is the other 
extreme of the range of galaxy activity compared to ellipticals; that 
is where we expect to find the most numerous and luminous X-ray 
binaries), but also a good sampling of all morphological types from 
early-type (large bulge) to late-type (small or no bulge) spirals and 
with a sufficiently large range of specific SFR (so as to differenti- 
ate between galaxies dominated by LMXBs, and those dominated 
by HMXBs). Each target should be observed deep enough to be 
complete at least down to 1039 erg s7! (in most cases, we required 
a detection limit of «3 x10°8 erg 871). Finally, the total exposure 
time has to be kept within plausible limits. 

There are about 70 Virgo spirals with SFR 20.3 Mo yr, but 
the precise list depends on the various alternative proxies for the 
SFR (Ha luminosity, or radio continuum, or far-IR, or far-UV plus 
24-um luminosity) as well as the adopted distance to the galaxies 
within Virgo. Moreover, some of the SFR studies we use for the 
analysis in this paper were not available at the time of our sample 
selection. With these caveats, we selected a representative sample? 
of about 60 galaxies with SFR 20.3 Мо yr"! , plus about 15 galaxies 
with a lower SFR, but which were either already in the Chandra 
archives or were in the same field of view of another target with 
a higher SFR. In total, we compiled a list of 75 spirals: a sample 
large enough to be representative of the population of actively star- 
forming galaxies in Virgo, and to provide a sizeable population of 
ULXs. In particular, we included all the «20 spirals with SFR >1 
Мо yr! . Of the 75 galaxies in the sample, 46 were observed for the 
first time by Chandra through our Large Program. For another six 
galaxies, new observations were combined with a previous snapshot 
Observation in order to reach the required detection limit, though 
any detected sources were subsequently analyzed separately for 
each observation. Finally, the remaining 23 galaxies already had 
sufficiently deep archival observations. In the rest of this Section 
we summarize the main properties of our sample. We will then 
describe our data analysis and main results, in particular for the 
ULX census. A more detailed analysis of the most interesting or 
peculiar individual sources will be presented in follow-up papers. 


22 Cluster membership 


Seventy-four out of our 75 galaxies (Figure 1) are included in the 
Virgo Cluster Catalog (VCC) of Binggeli et al. (1984), which is 
often used as the standard reference for a galaxy census in this 
cluster. The only galaxy in our sample that falls outside the VCC 
footprintis NGC 4713. However, itis included in the Extended Virgo 
Cluster Catalog (Kim et al. 2014), and it is a likely cluster member 
based on its recession speed and redshift-independent distance. 

A VCC classification does not of course guarantee that a galaxy 
is physically associated with Virgo. This is especially true for spirals, 
which are more spread out around the edges of the cluster, while 
ellipticals are more concentrated in the central region, around M 87 
and M 49. As a test that all 75 galaxies are at least plausible cluster 
members, we inspected their radial velocity distribution with respect 
to the cosmic background, using the values of vay reported in 


? We do not claim that the sample is complete within a certain threshold; 
simply that it is well representative of most spirals with high SFR. All the 
420 spirals with SFR z1 Мо уг”! are included. Among those with a SFR 
between ~0.4-1Mo yr-!, we are aware of the following six Virgo galax- 
ies missing from the sample: NGC 4294, NGC 4634, NGC 4409, IC 3476, 
NGC 4630, and NGC 4207. A case can be made also for NGC 4383 (with an 
uncertain classification between SO and Sa peculiar), and NGC 4651 (located 
at the northern edge of Virgo, probably a member). 
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Figure 3. Ratio between the median redshift-independent distance (from 
NED) and the conventional Virgo sub-structure distance (17, 23 or 32 Mpc), 
for the 75 galaxies in our sample. Three galaxies with very discrepant dis- 
tance measurements are labelled in the plot and discussed in Section 2.3. 
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Figure 4. Morphological type distribution of the galaxies in our sample, ex- 
pressed with the "Hubble stage" index T (horizontal axis) (de Vaucouleurs 
1959; de Vaucouleurs et al. 1991) and with the corresponding Hubble- 
Sandage-de Vaucouleurs classification (label over each histogram bin). The 
solid red histogram is for the galaxy morphologies reported in NED; the 
dash-dotted blue histogram is for the morphological classification of Hyper- 
LEDA. 


HyperLEDA?. We obtain a broad, structured, tri-modal distribution 
(Figure 2), which spans «3000 km 571, consistent with the findings 
of Binggeli et al. (1987) from the whole VCC. A few additional 
galaxies with Chandra archival data are physically located inside 
the VCC footprint (thus, they do have a VCC number) but their 
radial velocity is >3000 km s~!, which locates them most likely 
behind the cluster: we did not include those outliers in our sample 
of 75 candidate members. 


2.3 Distances 


Recession speed is not a useful proxy for the distance of individual 
Virgo galaxies, because it is comparable to the proper motion inside 
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Table 1. Main properties of the sample of 75 Virgo spirals 


Galaxy 


a) 


МСС 4178 
МСС 4192 
МСС 4197 
МСС 4206 
МСС 4212 
МСС 4216 
NGC 4222 
МСС 4237 
МСС 4254 
МСС 4276 
МСС 4293 
МСС 4298 
МСС 4299 
МСС 4302 
МСС 4303 
МСС 4307 
МСС 4312 
МСС 4313 
МСС 4316 
МСС 4321 
МСС 4330 
МСС 4343 
МСС 4356 
МСС 4380 
ІС 3322А 
МСС 4388 
МСС 4390 
ІС 3322 
NGC 4394 
МСС 4396 
МСС 4402 
МСС 4405 
NGC4411A 
NGC 4407 
NGC 4412 
NGC 4416 
NGC 4411B 
NGC 4419 
NGC 4424 
NGC 4430 
NGC 4429 
NGC 4438 
NGC 4445 
NGC 4450 
NGC 4451 
ІС 3392 
NGC 4457 
NGC 4469 
NGC 4470 
NGC 4480 
NGC 4492 
NGC 4496A 
NGC 4498 
NGC 4501 
NGC 4519 
NGC 4522 
NGC4527 
NGC 4532 
NGC 4535 
NGC 4536 
NGC 4548 
NGC 4567 
NGC 4568 
NGC 4569 
МСС 4571 
МСС 4579 
МСС 4580 
МСС 4606 
МСС 4607 
МСС 4639 
МСС 4647 
NGC 4654 
МСС 4689 
МСС 4698 
МСС 4713 


Other Names 


(2) 


VCC 0066 
M98, VCC 0092 
VCC 0120 
VCC 0145 
NGC 4208, VCC 0157 
VCC 0167 
VCC 0187 
VCC 0226 
M99, VCC 0307 
VCC 0393 
VCC 0460 
VCC 0483 
УСС 0491 
VCC 0497 
M61, VCC 0508 
VCC 0524 
VCC 0559 
VCC 0570 
VCC 0576 
M100, VCC 0596 
VCC 0630 
VCC 0656 
VCC 0713 
VCC 0792 
VCC 0827 
VCC 0836 
VCC 0849 
УСС 0851 
VCC 0857 
VCC 0865 
VCC 0873 
VCC 0874 
VCC 0905 
NGC 4413, VCC 0912 
VCC 0921 
VCC 0938 
VCC 0939 
VCC 0958 
VCC 0979 
VCC 1002 
VCC 1003 
VCC 1043 
VCC 1086 
VCC 1110 
VCC 1118 
VCC 1126 
VCC 1145 
VCC 1190 
NGC 4610, VCC 1205 
VCC 1290 
VCC 1330 
NGC 4505, VCC 1375 
VCC 1379 
M88, VCC 1401 
VCC 1508 
VCC 1516 
VCC 1540 
VCC 1554 
VCC 1555 
VCC 1562 
M91, VCC 1615 
VCC 1673 
VCC 1676 
M90, VCC 1690 
VCC 1696 
M58, VCC 1727 
VCC 1730 
VCC 1859 
VCC 1868 
VCC 1943 
VCC 1972 
VCC 1987 
VCC 2058 
VCC 2070 
EVCC 1189 


Type 


(3) 


SB(rs)dm 
SAB(s)ab 
Sd 
SA(s)be 
SAc 
SAB(s)b 

Sd 
SAB(rs)| 
SA(s)c 
Sc 
(R)SBO/a(s) 
SA(rs)c 
SAB(s)dm 
Sc 
SAB(rs) 
Sb 
SA(rs)ab 
SA(rs)al 
Scd 
SAB(s)bc 
Scd 
SA(rs)b 
Scd 
SA(rs)b 
SB(s)cd 
SA(s)b 
SAB(s)c 
SAB(s)cd 
(R)SB(r)b 
SAd 
Sb 
SAO/a 
SB(rs)c 
(R’)SB(rs)ab 
SB(r)b pec 
SB(rs)cd 
SAB(s)cd 
SB(s)a 
SB(s)a 
SB(rs)b 
50/а 
ЅА0/а(ѕ) pec 
Sab 
SA(s)ab 
Sbe 
SAb 
(R)SABO/a(s) 
SB(s)0/a 
Sa 
SAB(s)c 
SA(s)a 
SB(rs)m 
SAB(s)d 
SA(rs)b 
SB(rs)d 
SB(s)cd 
SAB(s)bc 
IBm 
SAB(s)c 
SAB(rs)bc 
SB(rs)b 
SA(rs)bc 
SA(rs)bc 
SAB(rs)ab 
SA(r)d 
SAB(rs)b 
SAB(rs)a pec 
SB(s)a 
SBb 
SAB(rs)bc 
SAB(rs)c 
SAB(rs)cd 
SA(rs)bc 
SA(s)ab 
SAB(rs)d 


6 


6 


R.A.(J2000) 


(4) 


12 12 46.45 
2 13 48.29 
2 14 38.55 
215 16.81 
2 15 39.36 
2155444 
21622.52 
21711.42 
2 18 49.60 
2 20 07.48 
221 12.89 
221 32.76 
2 21 40.55 
221 42.48 
221 54.90 
12 22 05.68 
2 22 31.36 
2 22 38.54 
12 22 42.24 
2 22 54.83 
12 23 17.25 
12 23 38.70 
12 24 14.53 
2252217 
22542.56 
22546.75 
22550.67 
122554.10 
12:25 55.53 
22558.82 
2 26 07.65 
22607.15 
2 26 30.10 
2 26 32.25 
12 26 36.08 
2 26 46.72 
12 26 47.23 
12 26 56.44 
12 27 11.61 
2272641 
2 27 26.51 
12 27 45.59 
12 28 15.93 
2 28 29.63 
12 28 40.55 
12 28 43.26 
12 28 59.0 
12 29 28.03 
2 29 37.78 
2 30 26.78 
12 3059.71 
12 31 39.2 
12 31 39.56 
12 31 59.16 
12 33 30.25 
12 33 39.7 
12 34 08.42 
12 34 19.33 
12 34 20.3 
12 34 27.05 
12 35 26.45 
12 36 32.7 
12 36 34.26 
12 36 49.79 
12 36 56.38 
12 37 43.52 
12 37 48.39 
12 40 57.54 
12 41 12.40 
12 42 52.39 
12 43 32.31 
12 43 56.58 
12 47 45.56 
12 48 22.91 
12 49 57.87 


Dec.(12000) 


ES 
ЕТ 


(5) 


05157.5 
454012 


+05 48 20.6 


+ 
+ 
+ 
+ 
+ 
+ 


3 01 26.3 
3 54 05.4 
3 08 57.8 
3 18 25.4 
5 1926.3 
424594 


+07 41 30.7 


4 
Ф 
+ 
$ 


8 22 56.6 
4 36 22.2 
1 29 59.7 
4 35 53.9 


+04 28 25.1 
+09 02 37.1 


ES 
+ 


5 3216.5 
148 03.4 


+09 19 56.9 


+ 


5 49 18.5 


+11 22 04.7 
+06 57 14.7 
+08 32 09.1 


+ 


0 01 00.5 


+07 13 00.0 


$ 
+ 


2 39 43.5 
0 27 32.6 


+07 33 17.2 
+18 12 50.6 


аб 
+ 
+ 


5 4017.3 
3 0648.0 
6 1051.6 


+08 52 20.0 


+ 


2 36 39.6 


+03 57 52.9 
+07 55 08.4 
+08 53 04.6 
+15 02 50.6 
+09 25 14.4 
+06 15 46.0 
+11 06 27.8 
+13 00 31.8 
+09 26 10.3 
+17 05 05.8 
+09 15 31.7 
+14 59 58.2 
+03 34 14.1 
+08 44 59.7 
+07 49 27.1 
+04 14 47.7 
+08 04 40.3 
+03 56 22.1 
+16 51 10.0 
+14 25 13.4 
+08 39 17.1 
+09 10 30.1 
+02 39 13.2 
+06 28 03.7 
+08 11 51.9 
+02 11 17.3 
+14 29 46.8 
+11 15 28.8 
+11 14 20.0 
+13 09 46.6 
+14 13 02.5 
+11 49 05.5 
+05 22 06.7 
+11 54 44.0 
+11 53 11.9 
+13 15 26.6 
+11 34 54.7 
+13 07 36.1 
+13 45 46.1 
+08 29 14.6 
+05 1841.1 


V3K 
(km 571) 
(6) 


714+2 
212+12 
2427+10 
036+2 
246+2 
468+3 
563+2 
192+1 
2737+8 
2959+1 
217+3 
456+4 
576+2 
1457+16 
912+2 
397+4 
475+3 
767+3 
1589+1 
89842 
1938+37 
133945 
432+6 
287+4 
335+3 
2845+3 
433+3 
1526+5 
1232+4 
208+6 
563+2 
2052+5 
1616+2 
426+4 
2603+18 
1729+2 
1608+2 
96+23 
776+3 
178842 
1325467 
42747 
684+3 
2273+3 
1195+2 
1992+2 
1228+2 
918+11 
2674+2 
2774+3 
2079+10 
2072+2 
1824+4 
2605+2 
1555+2 
2653+4 
2077+1 
2350+4 
2296+1 
2149+1 
802+2 
2588+8 
2563+6 
103+2 
657+3 
1842+2 
1372+2 
1969+7 
2580+8 
1302+7 
1733+2 
1362+4 
1924+5 
1333+2 
980+2 


dy 
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(Mpc) 


(8) 


13.20+3.00 


6.10+ 


25.30+4 
8.10+4 


3.00+ 


6.05+ 


3.55+2.93 
26.40+3.92 
8.70+3.40 
7.05+2.62 
5.50+2.08 
21.30+5.16 
20.40+5.19 
4.40+2.05 
27.90+* 

5.90+3.13 
5.80+2.54 
21.80+2.96 
7.40+8.60 
2.25+7.32 
21.95+3.39 
0.90--0.96 
4.15+2.92 
26.30+4.19 


M 


19.30--1.56 
26.80+4.02 
22.70+6.06 
9.3043.56 


12 
.96 


23.90+7.81 
20.15+2.58 
6.8043.13 


97 


4.50--3.83 
7.85+3.32 
6.8043.15 


40 


35.60+* 

42.10+* 
22.40+7.92 
16.10+3.62 
16.40+0.80 


5.85+2.77 
5.80+4.18 
11.30+3.54 
18.50+3.13 
5.25+3.36 
27.00+3.95 
13.65+3.86 
10.70+3.15 
6.75+0.07 
6.40+6.60 
42.20+3.80 


19.30+3.54 
15.70+1.14 
14.55+3.62 
17.10+4.58 
19.60--8.48 
16.05+4.73 
14.05+0.62 
13.00+3.45 
15.95+1.42 
14.95+1.73 
15.60+1.06 
22.45+5.79 
17.40+5.96 
11.45+3.16 
15.10+2.33 
16.80+3.07 
17.35+3.15 
14.10+4.31 
19.70+6.55 
21.95+1.97 
16.80+3.04 
14.40+3.15 
15.70+3.33 
22.15+6.77 
14.8043.55 


Іог(М. /Мо) 


(9) 


9.84+0.09 
0.63+0.09 
9.92+0.10 
0.02+0.10 
0.09--0.10 
1.00+0.09 
9.61+0.10 
0.25+0.10 
0.29+0.10 
9.81+0.10 
0.45+0.09 
0.05+0.09 
9.51+0.10 
0.40+0.09 
0.30+0.10 
0.38+0.09 
9.68+0.09 
0.01+0.09 
10.10--0.09 
0.71+0.10 
9.85+0.10 
10.33+0.09 
0.21+0.10 
0.41+0.10 
9.85+0.10 
0.10+0.09 
9.55+0.10 
9.36+0.10 
10.52+0.10 
9.12+0.10 
9.87+0.09 
9.86+0.10 
9.27+0.11 
9.79+0.10 
10.18+0.10 
10.29+0.10 
9.94+0.10 
10.31+0.09 
9.93+0.09 
9.66+0.10 
10.97+0.10 
10.49+0.10 
9.69+0.10 
10.75+0.10 
9.95+0.10 
9.75+0.09 
10.20--0.10 
10.42+0.10 
9.40+0.10 
10.42+0.10 
10.35+0.10 
9.62+0.10 
9.46+0.10 
11.03+0.09 
9.69+0.10 
9.75+0.10 
10.37+0.09 
9.28+0.10 
10.43+0.10 
10.22+0.09 
10.66--0.10 
10.28--0.10 
10.25+0.09 
10.53+0.09 
10.13+0.10 
10.96+0.10 
10.11+0.10 
9.77+0.10 
9.63+0.10 
10.29+0.10 
0 
1 
0 
0 
0 


10.010. 
10.120. 
10.36+0. 
10.87+0. 
9.46+0.1 


SFR(12um) 
(Мо угт!) 
(10) 


0.59+0.19 
2.03+0.66 
1.00+0.33 
0.42+0.15 
1.99+0.69 
1.66+0.56 
0.49+0.18 
1.30+0.47 
7.37+2.45 
0.54+0.20 
0.75+0.26 
1.39+0.47 
0.65+0.27 
1.58+0.55 
5.14+1.78 
0.64+0.22 
0.29+0.10 
0.33+0.1 
1.22+0.44 
7.45+2.54 
0.51+0.18 
1.02+0.37 
0.30+0.1 
0.6340.2 
1.10+0.39 
3.01+1.06 
0.42+0.15 
0.33+0.1 
0.54+0.19 
0.25+0.09 
1.33+0.44 
0.44+0.16 
0.14+0.05 
0.39+0.13 
2.34+0.8 
2.09+0.72 
0.37+0.13 
1.47+0.50 
0.55+0.18 
0.51+0.16 
0.16+0.06 
0.35+0.11 
0.20+0.07 
0.58--0.20 
1.01+0.36 
0.28+0.09 
0.41+0.13 
0.18+0.06 
0.42+0.14 
2.14+0.73 
0.24+0.09 
0.72+0.24 
0.35+0.12 
6.00--2.08 
.26+0.45 
0.43+0.15 
4.53+1.57 
0.60--0.19 
3.82+1.30 
3.64+1.29 
41+0.48 
2.18+0.75 
3.94+1.37 
.58+0.54 
:04+0.35 
.93+0.67 
0.48+0.17 
0.23+0.07 
0.96+0.34 
0.81--0.27 
:91+0.66 
2.93+1.01 
.34+0.46 
0.88+0.31 
0.73+0.24 


SFR(22um) 
(Mo yr!) 
(1) 


0.50+0.19 
1.36+0.51 
0.99+0.37 
0.28+0.11 
1.55+0.61 
0.93+0.36 
0.39+0.16 
0.93+0.38 
5.51+2.09 
0.35+0.15 
1.21+0.47 
0.92+0.36 
0.86+0.36 
0.96+0.38 
3.95+1.55 
0.48+0.19 
0.18+0.07 
0.20+0.08 
0.80--0.32 
6.05+2.35 
0.34+0.14 
0.69+0.28 
0.36+0.14 
0.35+0.13 
0.77+0.31 
5.11+2.03 
0.31+0.12 
0.19+0.07 
0.33+0.13 
0.18+0.07 
0.96+0.36 
0.30+0.12 
0.08::0.04 
0.27+0.11 
3.04+1.19 
1.34+0.53 
0.21--0.08 
2.05+0.79 
0.59+0.22 
0.32+0.12 
0.15+0.06 
0.24+0.09 
0.11+0.04 
0.40+0.16 
0.77+0.31 
0.18+0.07 
0.34+0.12 
0.11--0.04 
0.29+0.11 
1.68+0.65 
0.09+0.04 
0.62+0.21 
0.23+0.09 
4.19+1.65 
.33+0.54 
28+0.11 
92+1.53 
93+0.35 
34+1.30 
74+1.86 
88+0.34 
4340.56 
3.07+1.21 
.31+0.52 
0.65+0.25 
.29+0.50 
0.34+0.13 
0.14+0.05 
0 
0 


овкэсоос 


.73+0.29 
.58+0.22 
.36+0.53 
2.43+0.95 
0.86::0.33 
0.39--0.16 
0.63+0.24 


ЅЕЕ(В15) 
(Мо yr!) 
(12) 


0.51 
1.16 
МА 
0.28 
1.09 
0.66 
0.26 
0.48 
517 
0.33 
МА 
0.81 
NA 
0.86 
3.84 
0.22 
0.13 
0.11 
0.55 
5.46 
0.34 
0.40 
0.18 
0.27 
0.73 
3.01 
0.37 
0.15 
NA 
0.24 
0.65 
NA 
NA 
0.18 
240 
1.12 
NA 
1.21 
0.32 
0.25 
NA 
0.20 
0.11 
0.27 
0.63 
0.10 
0.29 
NA 
0.34 
1.42 
0.07 
0.80 
0.21 
3.98 
1.05 
0.34 
2.89 
1.18 
2.71 
3.49 
0.62 
0.65 
2.83 
0.89 
NA 
1.55 
0.20 
0.06 
0.62 
0.57 
1.07 
1.96 
0.54 
NA 
0.82 


Col. (3): from NED. Col. (6): radial velocity with respect to the cosmic microwave background (from HyperLeda). Col. (7): galaxy distance assuming the 
conventional structure of the Virgo Cluster, that is: 17 Mpc for the Virgo A subclump and the N, S and E clouds; 23 Mpc for the Virgo B subclump; 32 Mpc 
for the W and M clouds (as adopted, for example, by Gavazzi et al. 1999, Boselli et al. 2015, and the GOLDMine database). Col. (8): median 
redshift-independent distance from NED, and standard deviation of the published values. When a galaxy has Cepheid distance measurements, we used the 
median and standard deviation of those values only. Col. (9): stellar mass, derived from WISE W1 and W2 fluxes, as calibrated by Cluver et al. (2014) (see 
also Jarrett et al. 2019). The median NED distances were adopted. Col. (10): SFR from the WISE W3 flux, derived with the calibration from Cluver et al. 
(2017) (see Eq. 2), and assuming median NED distances. Col. (11): SFR from the WISE WA flux, derived with the calibration from Cluver et al. (2017) (see 
Eq. 3). Col (12): average SFR computed by Boselli et al. (2015), from four different proxies. We have rescaled their published values (originally based on 


sub-structure distances) to our median NED distances. The 1-c errors for the SFRs іп Col (12) аге «2096 of the values. 
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Figure 5. WISE classification of the galaxies in our sample, following 
Jarrett et al. (2019). All the galaxies follow the “star formation sequence” 
(dotted blue line) of normal galaxies, except for NGC 4388 (a Seyfert 2 
galaxy), in which the W/ colour is contaminated by the active nucleus. Spiral 
galaxies on the left of the sequence are more bulge-dominated, while those 
on the right have a younger and more active disk. The most bulge-dominated 
galaxy in our sample is NGC 4429 (classified as S0/a in HyperLEDA), while 
the most disky are NGC 4519 and NGC 4254 = M99. 
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Figure 6. Stellar mass distribution, from the mid-infrared luminosity (WISE 
colours), assuming the median redshift-independent distances from NED 
(Section 2.3). 


the cluster. There are two common alternative ways of allocating 
distances in Virgo cluster surveys. One method is to divide the clus- 
ter into sub-clusters and “clouds”, determine an average distance for 
each substructure, and attribute that distance to all the galaxies in the 
substructure. The most commonly used subdivision (Gavazzi et al. 
1999, 2002; Boselli et al. 2015) consists of a cluster A around M 87, 
located at an assumed (partly conventional) distance of 17 Mpc; a 
cluster B around the position of M 49, at 23 Mpc (but M 49 itself 
is only projected in front of cluster B, and is allocated a distance of 
17 Mpc); clouds N, S, E at 17 Mpc; clouds W, M at 32 Mpc (see 
Fig. 8 in Gavazzi et al. 1999 for a definition of the conventional sky 
boundaries of those regions). This is also the convention adopted 


in the Galaxy On-Line Database MIlano NEtwork (GOLDMine^, 
Gavazzi et al. 2003). 

Interms of our sample, 29 out of our 75 spirals belong to cluster 
A, 15 to cluster B, another 15 to cloud S, 8 to cloud N, 7 to cloud 
E and 1 to cloud W. Thus, in this classification, 59 of our galaxies 
have a distance of 17 Mpc, 15 are at 23 Mpc, and one (NGC 4197) is 
at 32 Mpc (see Table 1, Col. 7). The relative high number of targets 
in clouds S (in particular, at its southern end, for Declination £5?) 
and N is not surprising, as those two infalling subgroups consist of 
spirals for ~80% of their members (Gavazzi et al. 1999), while the 
ellipticals are more concentrated near the centre. 

While attractive for its simplicity, the main shortcoming of this 
classification is that the projected sky location inside a subgroup 
boundary is no guarantee that a particular galaxy is not in the 
foreground or background: the clouds and sub-clusters themselves 
may be several Mpc thick. An example of this issue is discussed in 
Boselli et al. (2002), who note the presence of gas-rich spirals in 
cluster A alongside the more common (and expected) gas-deficient 
ones. They suggest that the gas-rich galaxies are located «10 Mpc 
behind the core of Virgo, at the outskirts of the cluster, where 
they have not yet been ram-pressure stripped of their gas by the 
intracluster medium. One of the best examples of such gas-rich 
galaxies is NGC 4480. 

The alternative choice, adopted here, for the galaxy dis- 
tances is to use the database of redshift-independent measure- 
ments in the literature, as reported in the NASA/IPAC Extra- 
galactic Database (NED? ). Such distances come from a vari- 
ety of methods, and it is usually impossible to check the var- 
ious assumptions or limitations that affected the original anal- 
yses. When Cepheid distances were available (for 9 of the 75 
galaxies, namely NGC 4321, NGC 4424, NGC 4496A, NGC 4527, 
NGC 4535, NGC 4536, NGC 4548, NGC 4571 and NGC 4639), we 
used the median and standard deviation of those measurements. Oth- 
erwise, we adopted the median and standard deviation of all NED 
values (mostly Tully-Fisher measurements) after first inspecting the 
histogram of redshift-independent distances, removing outliers and 
re-deriving the median and its uncertainty. These are recorded in 
Table 1, Col. 8. Finally, we compared the median NED distances 
with the substructure distance (Figure 3). The scatter around a ratio 
of 1 is expected, given the typical standard deviation of several Mpc 
in the measured distances for each galaxy, and the intrinsic thick- 
ness of the substructures (also several Mpc). There are only three 
outliers (NGC 4412, NGC 4416, NGC 4480, labelled in Figure 3) 
which have Tully-Fisher distances (^40 Mpc) much larger than their 
substructure distances (17 Mpc). We do not have elements to decide 
at this stage which set of distances is more reliable, and whether 
those three galaxies are Virgo cluster members. The uncertainty in 
the distance to NGC 4480, in particular, affects our census of the 
most luminous ULXs, because the most luminous point-like source 
in that galaxy has either an X-ray luminosity ~ tor? erg 57! or barely 
above 1039 erg 87! depending on the choice of distance. 


2.4 Morphology and stellar masses 


Our sample includes a good balance of early- and late-type spirals 
(Figure 4). About half of the galaxies are of type SO/a to Sbc, and 
the other half from Sc to Im (with a little discrepancy between the 
classification in NED and HyperLEDA). There are eight Messier 


^ http://goldmine.mib.infn.it 
5 https://ned.ipac.caltech.edu 
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Figure 7. Solid red histogram: SFR distribution, inferred from the WISE 
12-um luminosity (Jarrett et al. 2019). Dash-dotted blue histogram: SFR 
distribution from the WISE 22-um luminosity with the scaling of Jarrett et 
al. (2019). 


galaxies (Messier 1781) in our sample: M58, M61, M88, M 90, 
M91, M98, M 99, M 100. 

An alternative classification of spiral galaxies is based on their 
mid-infrared colours, as measured by the Wide-field Infrared Survey 
Explorer (WISE; Wright et al. 2010). All 75 galaxies in our sample 
have been imaged and characterized with WISE: this is part of a 
larger project called the WISE Extended Source Catalog (WXSC). 
WISE photometry and global properties for the 100 largest galax- 
ies in the sky (first data release of the WXSC) was published by 
Jarrett et al. (2019): that sample already included a few of our 75 
Virgo targets. For all other galaxies, we used WISE results that will 
be published in follow-up WXSC data releases. For a detailed tech- 
nical description of how the total WISE luminosity of each galaxy 
was measured, see Jarrett et al. (2019). In short, the steps followed 
by the WSXC team were the following. First, Galactic foreground 
stars projected onto a galaxy were identified and removed. Then, 
a local background was computed from an elliptical annulus, with 
inner and outer radii located ~30% and ~50% beyond the maxi- 
mum extent of the galactic disk. Foreground stars were also masked 
from the background annulus. Each galaxy was assumed to have 
a constant ellipticity and to be axisymmetric. To estimate the total 
flux of the galaxies, Jarrett et al. (2019) simultaneously fit two Sér- 
sic functions (bulge + disk) to the surface brightness profile. See 
Jarrett et al. (2019) for details. 

WISE has four channels: W1, W2, W3 and WA, correspond- 
ing to the 3.4-шп, 4.6-um, 12-шп and 22-um bands, respectively. 
Plotting the W1 - W2 colours against the W2 — W3 colours defines 
a "star formation sequence" of normal galaxies (Jarrett et al. 2017, 
2019). All our galaxies are located along that sequence (Figure 5), 
except for NGC 4388, which is «0.3 mag above the sequence, in 
a region typically occupied by Seyfert and low-luminosity AGN 
(Mingo et al. 2016; Huang et al. 2017), which is indeed the case for 
NGC 4388. The reddest galaxy along the mid-infrared star forma- 
tion sequence is NGC 4429 (classified as S0/a in HyperLEDA); the 
bluest one is NGC 4254 = M 99 (an Sc galaxy). 

Stellar masses (M+) were derived from the mid-infrared lumi- 
nosity. Specifically, we used the relation from Cluver et al. (2014), 


6 https://vislab.idia.ac.za/research-wxsc 
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Figure 8. Comparison between the SFRs derived from the 12-um lumi- 
nosities, the 22-j2m luminosities, and the average SFRs from Boselli et 
al. (2015), B15. The dashed line is the one-to-one relation. On average, the 
22-um SFRs are ~20% higher than the B15 values (plotted in blue) and 
2:204 lower than the 12-р values (plotted in red). The difference between 
the three estimates is relatively more prominent for low SFRs. Typical uncer- 
tainties on the WISE SFRs from Equations 2 and 3 are around 33-38 percent 
(see Table 1), and around 20 percent for the B15 values. The red datapoint 
in the top left sector of the plot and the blue datapoint in the bottom right 
sector are labels, displayed to represent the typical size of the error bars. 


based on the WISE photometry; see also Jarrett et al. (2013, 2019): 
log (М, / Іду) = -2.54(W1 - W2) – 0.17, (1) 


where M. is in solar units, W1 апа W2 аге the magnitudes in the 
3.4-шп and 4.6-шп bands, respectively, and Lw | is the luminosity 
in the W1 band in units of solar luminosity. (The absolute magnitude 
of the Sun in the W1 band is My ,o = 3.24 mag). This definition 
of M. is in excellent agreement (better than 10%) with the mass 
relation used by Graham et al. (2019) in their study of Virgo spirals, 
namely М. /Г.к = 0.6, where Lx is the K-band luminosity. (The 
good agreement between the K-band and WISE mass values holds 
also for NGC 4388, the galaxy where the WISE photometry was 
moderately affected by an AGN, as noted earlier). The main reason 
we adopt the WISE scaling rather than the K-band scaling in this 
paper is because we have also used the WISE band luminosities as 
a proxy for the SFR (Section 2.5). 

From the mass values listed in Table 1 (Column 9), we see that 
~3/4 of our sample galaxies have masses between «5 x10? Мо 
and «5 x10!9°Mo (Figure 6). The four most massive galaxies 
are NGC 4216, NGC 4429, NGC 4501 (= M88) and NGC 4579 (= 
М 58): all four of them have М, ~ 1 x 10!! Мо. The mass values 
scale with the square of the assumed distances. The values listed 
in Table 1 and adopted for subsequent analysis are for the median 
NED distances, after correcting for outliers. The error range for each 
galaxy includes only the photometric error in WISE W1 and W1 — 
W2, not the systematic uncertainty in the distance; it also does not 
include the scatter around the best-fitting relation, which we estimate 
as +0.15 dex, from Fig. 6 of Cluver et al. (2014). In total, the stellar 
mass of our galaxy sample is М. tot = 1.5 х 1012 Mo, if we assume 
the redshift-independent NED distances, ог М, tot = 1.7 х 10!2M о 
if we use the substructure distances. The main reason for computing 
a total stellar mass is that we shall estimate the number of LMXB 
ULXs per unit stellar mass in our sample of star-forming galaxies 
(Section 4.4, Section 5), and compare it both with the value found in 
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the early-type galaxies of the AMUSE sample (Plotkin et al. 2014) 
and with theoretical predictions of X-ray binary populations. 


2.5 Star formation rates 


To characterize the star formation properties of the sample, we 
looked for SFR proxies that were available for all 75 galaxies, for 
consistency. The luminosity in the WISE bands W3 (12 um) and 
WA (22 um) can be used as a proxy for the total infrared luminosity 
and, therefore, for the SFR (Jarrett et al. 2013; Cluver et al. 2014, 
2017; Jarrett et al. 2019). The relations calibrated by Cluver et al. 
(2017) are: 


SFR 
log | ——— 
Moyr! 

SFR 
log | ——— 
Moyr-l 


for the W3 and W4 bands, respectively. In Equations (2-3), v is the 
central frequency of the band, and the luminosity vL, is normalized 
to the bolometric luminosity of the Sun (Lo = 3.83 x 1033 erg 
51), not its band-limited luminosity. The contribution of the stellar 
continuum was subtracted from Lj24m and Гәзит as described 
in Cluver et al. (2017). A Kroupa initial mass function (IMF) was 
adopted (Kroupa 2001). 

The WISE SFRs are presented in Table 1 (columns 10 and 11) 
and are graphically shown in the upper panel of Figure 9. The total 
SFR урт of the sample is (104 + 6)Mo yr !, if we adopt redshift- 
independent distances, or (115 + 8)Mo уг”! for substructure dis- 
tances. The total ЗЕ т is (86 +6)Mo уг”! or (95--7)Мо yr! 
in the two cases. (Note that the two sets of WISE-derived SFRs scale 
as 41-75 and d!-89, respectively, as a function of assumed galaxy 
distances, while stellar masses scale as а), It is clear (Figure 7) that 
there is a systematic offset in the distribution of SFR }2,,m compared 
with the corresponding values of БЕЕ2 т. This is evident in par- 
ticular at lower SFRs (x0.3 Мо уг” !), where SFR рт tends to be 
7.1.52 times higher than SFR29,,m; at higher SFRs (22 Mo уг), 
the excess of SFR}2,m over SFR224m is only ~10%. Cluver et al. 
(2017) and Jarrett et al. (2019) argue that, for normal galaxies (ex- 
cluding extreme dust-obscured starbursts, which are not a concern 
in our sample), SFR,2,, is a better proxy of total infrared lumi- 
nosity and SFR than SFR22,,m, because it is less contaminated by 
polycyclic aromatic hydrocarbon emission and less dependent on 
metallicity. 


Il 


vL 
(0.873 + 0.021) | =] (2) 
(9) 


2(7.62 + 0.18) 


Il 


vL 
(0.900 + 0.027) | ae) (3) 
о 


-(7.87 + 0.24) 


In addition to the far-infrared luminosity, the most common 
empirical proxies for the SFR in the literature are the Ha luminosity 
(corrected for extinction), the far-UV plus 24-um luminosity, and the 
20-cm radio continuum luminosity; see Kennicutt (1998) for a well- 
known review. Boselli et al. (2015) (B15) collected a large database 
of SFR measurements in nearby spirals, determined with those three 
methods, as part of their Herschel Reference Survey; some of the 
SFR values in B15 came from their own Ha observations and data 
analysis, others were taken from published literature. When more 
than one proxy was available for a galaxy, B15 also included the 
average between the different values. Their SFR catalog includes 64 
of our 75 galaxies (Table 1, Column 12). 

The list of average SFRs in Boselli et al. (2015) does not in- 
clude errors. We estimated an approximate error from the dispersion 
of the SFR values reported for the individual proxies around their 
respective average SFRs, for all galaxies with more than one SFR 
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Figure 9. Top panel: 12-um SFR versus M... Loci of constant sSFR are 
marked with dashed lines. Bottom panel: as in the top panel, for the average 
B15 SFR. 


proxy. The SFRs from the different proxies have an approximately 
Gaussian distribution around the (normalized) average values, with 
a standard deviation с of about 20% of the mean. Thus, we take 
this standard deviation as an estimate of the 1-с error for each B15 
average SFR value. 

We compared the B15 average values with the WISE values 
(scaled to the same redshift-independent distance). The B15 values 
are lower than both the 12-~m and 22-um values (Figure 8). For 
the 64 galaxies that are in common between the WISE catalog 
and the B15 catalog, we determine a total 5ЕКв|5 = 66Mo уг, 
compared with SFR224m ~ 81Mo yr! and SFRj2, > 98Мо 
yr-!. The discrepancy occurs at all values of SFRs, although it 
is more accentuated at the low end (Figure 7). With the plausible 
assumption that the same fractional difference would be found in the 
remaining 11 galaxies, we conclude that the three SFR proxies used 
by B15 suggest a total SFR ~ 70Мо уг”! for our full sample. The 
large difference cannot be explained with the different choice of IMF 
(a Salpeter IMF in B15, Salpeter 1955); a conversion to a Kroupa 
IMF in B15 would reduce rather than increase the amount of stellar 
mass required to produce the same ionizing flux. An investigation 
into the different calibrations of the various SFR proxies is beyond 
the scope of this paper. 

Knowing the mass and the SFR, we can then determine the 
distribution of our sample galaxies in terms of their sSFR. This is 
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Figure 10. Solid red histogram: sSFR distribution, inferred from the WISE 
12-ит luminosity. Dash-dotted blue histogram: sSFR distribution from the 
WISE 22-um luminosity. 


a parameter that will be particularly useful later, when we interpret 
the X-ray luminosity of the galaxies. In the mass versus SFR] 2m 
plane (Figures 9,10), about half of our galaxies have sSFR < —10, 
which characterizes the “blue cloud” of star-forming galaxies (com- 
pare with Fig. 2 in Lehmer et al. 2019), and the other half can be 
described as “green valley” galaxies. If we adopt the B15 SFR, 
about two thirds of the galaxies are in the green valley. The sSFR 
distribution inferred from БЕЕ ит is of course intermediate be- 
tween the other two. From this rough classification, we can already 
speculate that on average, we should find comparable contributions 
to the X-ray luminosity from galaxies dominated by young HMXBs 
and those dominated by intermediate-age or old LMXBs (see Fig. 6 
in Lehmer et al. 2019 for the predicted transition between the two 
classes). We investigate this point with our Chandra study. 


3 X-RAY DATA ANALYSIS 


We selected and downloaded from the Chandra archives all the 
Advanced CCD Imaging Spectrometer (ACIS) observations that 
included our sample galaxies. That is a total of 110 observations 
from 2000 February to 2020 March: 106 observations placed the 
target on the back-illuminated ACIS-S3 chip, and the other 4 on 
the front-illuminated ACIS-I chips. Fifty-two of the observations 
were specifically obtained for this project (Proposal 18620568, PI 
R. Soria, observing cycle 18), with exposure times between 8 and 
15 ks, for a total exposure time of 551 ks; they all placed the target 
on S3. The other 58 were gathered from the archives, with various 
durations, from as short as 1.4 ks to as long as 148 ks, adding 
another 1,413 ks (Table 2, Columns 7-10). 

We reprocessed and analysed each dataset with the Chandra 
Interactive Analysis of Observations (сіло) software version 4.12 
(Fruscione et al. 2006), with calibration database version 4.9.1. In 
particular, we created new level-2 event files with the ciao tasks 
chandra repro. For galaxies with multiple observations, we applied 
reproject obs to align and create stacked images. We created indi- 
vidual and stacked images in the soft band (0.3-1.0 keV), medium 
band (1.0-2.0 keV), hard band (2.0-7.0) and full band (0.3—7.0 keV). 
We used the ps9 imaging tool (Joye & Mandel 2003) for displaying 
images, defining source and background regions, and performing 
various other tasks (e.g., smoothing, creating contour plots, align- 
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ing and overplotting X-ray contours over optical images, finding 
centroids of point-like sources, etc.). 


3.1 Integrated galaxy luminosities 


First, we measured or constrained the total X-ray luminosity of the 
X-ray binary population of their host galaxies (inside their D25), 
including the contribution from faint, unresolved sources. The main 
hurdles to overcome for such analysis were the choice of suitable 
background regions (especially in cases when the D»s of a galaxy 
occupied most of the S3 chip) and the potential additional contri- 
bution to the X-ray emission from diffuse hot gas and/or from an 
active nucleus. 

We have taken the B-band 25 mag arcsec ? isophotal radii, 
and thus diameters (025), from de Vaucouleurs et al. (1991). We 
drew corresponding elliptical regions with ps9, in a format suitable 
to сіло analysis. We then ran the сіло task specextract to create 
spectra and associated response and ancillary response files. More 
specifically, because we were extracting spectra of extended regions, 
we built spatially weighted ancillary response functions (specextract 
parameter “weight = yes") without a point-source aperture correc- 
tion ("correctpsf = no"). For each target galaxy, background regions 
on the same CCD were chosen as large as possible, outside the D25 
of the galaxy, but we avoided regions of low transmission (caused by 
molecular contamination) near the edges of the CCD array. Each of 
the six ACIS-S chips is 8’ x 8’, providing enough galaxy-free field- 
of-view. As noted previously, the bulk of the galaxy sample had 
their centre positioned on the S3 chip. We also built response and 
ancillary response files for the background regions with specextract. 

We adopted two methods to subtract the background and model 
the net galaxy emission. The first method consists of first modelling 
the background spectrum, scaling that model to the area of the 
galaxy, then holding this background model component fixed and 
adding one or more spectral components representing the emission 
from the galaxy. The second method simply subtracts the back- 
ground spectrum (suitably scaled) directly from the galaxy-plus- 
background spectrum. 

In the first method, the background model consists of a celestial 
X-ray background (which includes unresolved cosmic background 
from distant AGN, and possible intracluster emission from diffuse 
hot gas in Virgo) and a particle-induced instrumental background. 
The celestial X-ray background is modelled as a two-temperature 
thermal plasma plus a power law with spectral parameter values 
given by the Athena Wide Field Imager (Meidinger et al. 2020) 
background preparation document’, provided by Arne Rau (Max 
Planck Institute for Extraterrestrial Physics). We scaled the normal- 
ization of this component by the area of the background region on 
the sky. The instrumental background was modelled as a blackbody 
and three power-law continua superposed with three Gaussians. 
All parameters for the black body and the power-law components 
were allowed to vary in the fitting; however, the Gaussian energies 
and widths were fixed following the analytical particle background 
model of Bartalucci et al. (2014), and only their normalizations 
were allowed to vary. The best-fitting background model was then 
applied to the source spectrum, with all its parameters frozen, and 
a global normalization factor scaled by the ratio of the source area 
to the background area. 

Models were fitted to the spectra with the xsPEc (Arnaud 1986) 


7 https://www.mpe.mpg.de/ATHENA-WFI/public/resources/background/WFI- 
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package, version 12.11.0. The Cash statistics (Cash 1979) was used 
as the fit statistics (cstat in xsPEC). An absorbed power-law com- 
ponent was introduced to represent the galactic emission from (re- 
solved and unresolved) point-like sources. Initially, the absorption 
column density Ny of the galactic power-law component was fixed 
to the Milky Way line-of-sight value, and the power-law photon in- 
dex was fixed to Г = 1.7. In the event of an unacceptable fit, Мн was 
allowed to vary first, and then, if necessary, the power-law photon 
index was also thawed. 

The second method was a standard, direct subtraction of the 
background photon spectrum from the source spectrum (і.е., with 
the source and background spectra loaded in хѕрес as data file 
and bkg file, respectively). In this case, we did not need to assume 
any model for the background. We simply modelled the residual 
emission with an absorbed power-law, with free Ny and I' when 
we had 2 100 net counts, or a column density fixed to the Galactic 
line-of-sight value and Г = 1.7 for galaxies with very few net counts. 

We compared the results of the two methods, and found that 
they were consistent within the uncertainties. In this paper, we have 
reported the results from the first method. All fits spanned the 0.5- 
6 keV energy range, because the detected photons are strongly domi- 
nated by background above those energies. Observed (model) fluxes 
and de-absorbed emitted luminosities were then extrapolated and 
quoted over the 0.5-8 keV range, for consistency with the range 
chosen by Lehmer et al. (2019), which we will use extensively in 
comparison to our results. 


3.1.1 Caveat I: galaxies larger than the S3 chip 


In three cases, a sector at the outskirts of ће D25 region was 
outside the detector's field of view (namely, about half of the О25 of 
NGC 4178 and NGC 4579, and about 10% of NGC 4501). In those 
cases, the net galaxy flux and luminosity were scaled under the 
assumption of a uniform surface brightness. The resulting flux may 
be overestimated by this assumption, because the surface brightness 
typically decreases towards larger galactic radii. However, given the 
small number of galaxies affected by this issue, this likely over- 
estimation in flux does not significantly affect our population results; 
in fact, the uncertainty caused by the rescaling is within the overall 
model flux uncertainties for each of those galaxies. In another four 
galaxies with a large apparent size, the central part of the D»s 
falls onto the ACIS-S3 chip but a small sector falls onto either 
of the two adjoining chips (S2 or S4). In those cases, we defined 
local background regions that also spanned two chips, with the same 
fractional area as the source regions. We compared this method with 
the alternative possibility of computing the galaxy emission only 
over the fraction of the D5 inside the S3 chip, and then rescaling 
the total flux as explained above. The two approaches give the same 
result within the model uncertainties. 


3.1.2 Caveat 2: galaxies with a nuclear X-ray source 


A few of our galaxies contain a bright AGN, with 0.5-8 keV lu- 
minosity ~10°9-10*! erg 871: they are classified in the litera- 
ture as Seyferts (NGC 4388, NGC 4569, NGC 4579, NGC 4639, 
NGC 4698) or LINERs (NGC 4438, NGC 4450, NGC 4548). Those 
strong nuclear X-ray sources were immediately excluded from the 
extraction of ће D5 spectra of their host galaxies, with the use of 
"exclude" circular regions around the nucleus. In few cases, bright 
AGN readout streaks were visible on a chip and were also excluded 
from our source extraction. In all other cases, when we detected 


candidate nuclear sources but with an X-ray luminosity 510% еге 
571, we included the nuclear region іп ће 25 source extraction 
and subsequent model fitting. The emission of those faint, candidate 
nuclear point sources was extracted and modelled separately, and 
subtracted from the total D»5s5 luminosity of their host galaxies a 
posteriori. This is because modelling a weak nuclear source on its 
own and then subtracting its luminosity gives a more accurate result 
than simply placing an exclusion circle around the nuclear position. 
In total, 39 of the 75 galaxies (including the 8 mentioned above) 
have a point-like Chandra source within ~1’’.5 of the position of 
the optical nucleus reported in the literature. A separate paper on the 
nuclear X-ray properties of those galaxies suspected of hosting an 
intermediate-mass black hole can be found in Graham et al. (2021), 
and we have already presented preliminary results for candidate 
intermediate-mass nuclear BHs (based on archival Chandra data 
prior to our Large Program observations) in Graham et al. (2019). 


3.1.3 Caveat 3: galaxies with strong thermal-plasma emission 


In a few cases, a simple power law model was statistically unaccept- 
able. Inspection of the fit residuals (and, sometimes, direct inspec- 
tion of the images) indicated that a soft thermal component was also 
present іп the D25 emission, with typical emission features around 1 
keV, 1.3 keV and 1.8 keV; not surprisingly, this happened in galaxies 
with a high SFR (a beautiful example is NGC 4303). For those galax- 
ies, an additional apec model component was added. The flux from 
this thermal component was excluded from the galaxy flux reported 
here, because for this work we are specifically looking at the X-ray 
binary contribution (both the individually detected sources and the 
faint, unresolved ones). We are aware that part of the unresolved 
emission from a galaxy can be in the form of thermal plasma but 
come from point-like sources; for example, from accretion columns 
in magnetic Cataclysmic Variables (e.g., Revnivtsev et al. 2007; 
Li & Wang 2007; Li et al. 2011; Hong et al. 2016). However, such 
contribution is negligible compared with the power-law compo- 
nent produced by the X-ray binary population; it is also typically 
characterized by a temperature 210 keV, so that its spectrum is in- 
distinguishable from a hard power-law component, given the small 
number of counts for most of our galaxies, and the low Chandra 
sensitivity above 5 keV. It is only the thermal plasma component at 
temperatures ~0.3—1 keV, mostly from gas shock-heated by young 
SNRs, that (occasionally) stands out in a galaxy spectrum and can 
be separated from the power-law component. We leave a discussion 
of the hot gas emission to follow-up work. 


3.1.4 Caveat 4: background AGN projected onto the D25 


Area-scaling for background subtraction assumes that the cosmic 
background emission is uniform between source and background re- 
gions. This applies in principle both to unresolved emission and (on 
a Statistical level) to point-like sources (background AGN), leaving 
aside the issue of cosmic variance. This may not be true for indi- 
vidual galaxies, if their D25 ellipse includes опе or few bright AGN 
projected through the galaxy, which may not be exactly balanced by 
a scaled number of similar AGN in the background region, given 
the small-number statistics. We will illustrate the expected relative 
contribution of resolved background AGN in Section 4.2, when we 
discuss the ULX LF. Here we can anticipate that such a contamina- 
tion is not a significant problem. We will show that most galaxies 
have 0.5-8 keV X-ray luminosities z 1099, and we expect to find a 
total of only ~10 AGN projected inside the D25 of the 75 sample 
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galaxies, at a flux greater than this apparent luminosity (using the 
median NED distance for each galaxy). 


3.2 Point-source X-ray luminosities 


For this part of the analysis, we used again standard сіло tasks. 
In particular, we used wavdetect to identify point-like sources. We 
then checked and improved the centroid position with dmstat. In 
a few cases, especially when a source was observed far from the 
aimpoint, we further refined its position using a 2-D Gaussian fit 
with the sHERPA package (Freeman et al. 2001). For regions with 
significant diffuse emission from hot gas, we used the hard band 
or the 1.5—7 keV band to filter out the thermal plasma emission 
and identify accreting point-like sources. This includes also the 
candidate nuclear sources, which will be discussed in a follow-up 
paper. In this work, we present the results for the off-nuclear ULX 
population. 

We measured absorbed and absorption-corrected fluxes of ev- 
ery moderately bright point-like source (1.e., every source with an 
observed flux z 10-17 erg cm ? sl, corresponding to a luminosity 
of about a few x 1098 erg 5—1) with srcflux, over ће 0.5-8 keV band. 
In doing so, we used the ancillary response function to determine 
the shape of the point spread function (PSF) at the location of the 
source (і.е., “PSF Method = arfcorr" in srcflux). As usual for Chan- 
dra/ACIS studies, the size of the source extraction region depended 
on how far from the aimpoint each particular source was. For iso- 
lated sources near the aimpoint, the typical source extraction radius 
was 2" .5, or 2".0 for sources in more crowded regions. For back- 
ground subtraction, we defined local background regions at least 
four times larger that the corresponding source regions. 

We made two assumptions in our initial srcflux analysis of 
the ULX population: a) that the spectrum was a power law with 
a photon index Г = 1.8; b) that the absorbing column density 
Мн was limited to the line-of-sight Galactic value, obtained from 
HI4PI Collaboration et al. (2016) via NASA's High Energy Astro- 
physics Science Archive Research Center (HEASARC). The ratio- 
nale for the first assumption is that Г ~ 1.8 is the mean photon index 
over the 0.5—8 keV band determined from the Swartz et al. (2004)'s 
Chandra survey of ULXs in nearby galaxies. It is also the slope used 
for the conversion from count rates to luminosities in the follow-up 
work of Swartz et al. (2011). The effect of different choices of Г is 
small, over a plausible range of values. For example, for Cycle-18 
observations, with Ny = 3 x 1020 cm 2, the choice of Г = 1.6 
would increase the inferred unabsorbed luminosity by about 6%; 
the choice of T = 2.0 would decrease it by about 3%. The reason for 
the second assumption (line-of-sight Galactic absorption) is that we 
wanted to start from a ULX list as simple and model-independent 
as possible. 

In the second step of our point-source analysis (Section 4.2.1), 
we selected all sources with 250 counts. We used specextract to 
create spectra and associated response and ancillary response files 
for those sources. We then rebinned the spectra with the rroors 
(Blackburn 1995) task grppha, to 1 count per bin, and modelled 
them with xsPEc (Arnaud 1986), using the Cash statistics (Cash 
1979). We assumed a power-law model whenever possible, but un- 
like for our initial srcflux estimate, this time we left the photon index 
and the intrinsic Мн as free parameters. We used the cflux convolu- 
tion model in хѕрЕС to measure absorbed and de-absorbed fluxes in 
the 0.5-8 keV and 0.3-10 keV bands. 
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Figure 11. Luminosity of the X-ray binary population (0.5-8 keV band) 
of each galaxy, normalized to the SFR of the galaxy, plotted versus the 
specific SFR. Error bars are omitted for clarity, except for three galaxies 
in each plot, which visualize representative values of the error bars for 
low, intermediate and high sSFRs. Errors for all individual galaxies can be 
derived from the values of M., SFR and Lx listed in Table 1, Columns 
9-12, and Table 2, Column 2. Top panel: SFR values from the WISE 12- 
ит band; bottom panel: SFR values from B15 (hence, not all 75 galaxies 
have a datapoint in the bottom panel). The dashed diagonal red line is the 
theoretical luminosity expected for a pure LMXB population (і.е., scaling 
with the stellar mass), in the limit of perfect sampling of the LF (Lehmer et al. 
2019). The dashed horizontal blue line is the theoretical luminosity for a 
pure HMXB population, again in the limit of perfect sampling (Lehmer et al. 
2019). The red star marks the total value of Lx/SFR versus sSFR for all the 
sample galaxies. The red square is the expected value of Lx /SFR versus sSFR 
for the whole sample (i.e., the sum of the LMXB and HMXB luminosities 
predicted by Lehmer et al. 2019 at that value of sSFR). The green dash- 
dotted curves are the predicted median values of the observed luminosity, 
for a population of galaxies of various SFRs and stellar masses, based on 
the Monte Carlo simulations of Lehmer et al. (2019); the number next to 
each curve represents the stellar mass along that track. The green curves 
differ from a simple sum of LMXB plus HMXB luminosity, because of 
the incomplete sampling of the LF in each galaxy, especially those of lower 
masses. A few outliers with an exceptionally high Lx in relation to their SFR 
and stellar mass are labelled in both panels and will be discussed individually 
in follow-up work. 
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4 MAIN RESULTS 


We summarize our main results for the de-absorbed X-ray luminos- 
ity of all the sample galaxies, and the fluxes and luminosities of 
their ULXs, in Table 2. 


4.1 Galaxy luminosities 


The galaxy luminosities (Table 2, Column 2) include the ULXs 
but exclude the thermal plasma emission and the nuclear sources. 
This is because here we want to study the relation between the 
galaxy properties (stellar mass and SFR) and their X-ray binary 
populations. For galaxies with multiple observations of comparable 
duration, the luminosity values provided in Column 2 are an average 
luminosity. For galaxies with multiple observations of very uneven 
duration, only the deepest observations were used. When a galaxy 
contains ULXs, the ULX luminosity itself is a significant fraction of 
the total luminosity (/.e., the total point-source luminosity is usually 
dominated by the sources at the upper end of the LF). 

We stress that for this paper, our main interest is in the pop- 
ulation properties of the sample rather than in individual galaxies: 
specifically, whether the Virgo spirals are dominated by LMXBs 
or HMXBs?. The combined X-ray luminosity of LMXBs is pro- 
portional to the stellar mass (Gilfanov 2004; Lehmer et al. 2010; 
Boroson et al. 2011; Zhang et al. 2012; Lehmer et al. 2019), al- 
though other factors such as population age and specific frequency of 
globular clusters also affect the correlation (Zhang et al. 2012). The 
combined HMXB luminosity is proportional to SFR (Grimm et al. 
2003; Ranalli et al. 2003; Lehmer et al. 2010; Mineo et al. 2012a; 
Lehmer et al. 2019). Regardless of the precise value of the nor- 
malization constants, we expect a total X-ray luminosity Lx = 
а мхв M. + Внмхв SFR (Lehmer et al. 2010, 2019). Hence, if 
a population of galaxies are dominated by LMXBs, we expect 
Lx/SFR œ M,/SFR = 1/sSFR, that is we expect to find them 
along a line of slope —1 in the log (Lx /SFR) versus log sSFR 
plane. Instead, if the galaxies are dominated by HMXBs, we expect 
them scattered along a horizontal line, because Lx /SFR « constant. 
In the analysis of Lehmer et al. (2019), the transition from LMXB 
to HMXB domination occurs at sSFR ~ 107195 yr. Their best- 
fitting relation (which we will use for comparison with our data) is 


Ly- |52520 x SFR + 1.803 x М..о) x10? ergs}, (4) 


where the SFR is in units of Мо уг”! and M. 10 is the stellar mass 
in units of 1010 Мо. We have ignored here the small dependence of 
the SFR coefficient on the metal abundance (Lehmer et al. 2021); 
we leave this discussion to further work. 

Our analysis shows (Figure 11) that there is indeed a (1/sSFR) 
trend in the distribution of our galaxies at the lower end of the sSFR 
scale (where we expect LMXBs to dominate). This is more evident 
when we use the B15 values of SFR (Figure 11, bottom panel): in 
this case, about 1/3 of the 75 galaxies are in the region ofthe diagram 
in which the X-ray luminosity is dominated by LMXBs. This does 
not mean that 1/3 of the total X-ray luminosity of the Virgo spirals 
comes from LMXBs, because the galaxies dominated by HMXBs 
are also the most luminous and the most likely to contain ULXs. 

For sSFR > 10-105 yr, the galaxy distribution is very scat- 
tered and flatter than the (1/sSFR) trend, but most galaxies fall 


8 The possibility of off-centre, wandering intermediate-mass BHs among 
the point-source population will be discussed in further work. 


below the median tracks predicted by the Monte Carlo simula- 
tions of Lehmer et al. (2019) (despite a few exceptions with higher- 
than-expected X-ray luminosity, namely NGC 4532, NGC 4396A, 
NGC 4396 and NGC 4197). We had previously shown (Figure 6) 
that about two thirds of our sample galaxies have a stellar mass 
in the range 9.5 € log (М./Мо) x 10.5. However, the median 
point-source luminosity of our sample of galaxies in the high 
sSFR region (Figure 11) falls below the simulated median line for 
log (M«./ Мо) = —9.5 in Lehmer et al. (2019). Thus, the X-ray lumi- 
nosity of several of the most actively star-forming disks appears un- 
derestimated, compared with the scaling relations of Lehmer et al. 
(2019). Those scaling relations were derived from galaxies with 
discs orientated close to face-on, with low or intermediate incli- 
nation; however, our sample contains several galaxies seen at high 
inclination, with much higher column densities through their discs 
(particularly in the regions where HMXBs are preferentially lo- 
cated). 

On the same log (Lx/SFR) versus log sSFR plane we can 
also plot the total values from the whole sample, knowing that 
Lx tot & 3.6 х 104! erg 871 (from a simple addition of the 0.5-8 
keV luminosities of each galaxy) and the total SFR ~70-104 Mo 
yr-!. The total observed value of Lx /SFR (marked by a red star in 
the two panels of Figure 11) is lower than predicted by the model of 
Lehmer et al. (2019) (red squares in the two panels of Figure 11), 
when we include both the predicted HMXB contribution (dashed 
blue line) and the predicted LMXB contribution (dashed red line). 
In particular, the observed luminosity is a factor of 2 lower than 
predicted if we adopt the 12 set of SFRs (Figure 11, top panel), 
or a factor of 1.6 lower if we use the B15 set of SFRs (Figure 11, 
bottom panel). 

At this stage of our analysis, there is no glaring single reason 
for this discrepancy between inferred and predicted luminosities. 
We have already hinted at an underestimation of the de-absorbed 
luminosity from galaxies seen at higher inclination as the most 
likely cause. We will investigate this possibility further when we fit 
the individual spectra of the most luminous ULXs and model their 
luminosity distribution (Sections 4.2.1, 4.3). We will show there that 
the median absorption column density of the most luminous ULXs 
is Ny ~ 3x10?! спі 2. If we assume this column density, instead of 
line-of-sight Galactic absorption, our inferred luminosities would be 
increased by a factor of ~1.5 (for observations in Cycle 4 or earlier), 
a factor «1.4 (observations in Cycles 5-14) and a smaller factor for 
more recent observations (~1.2 for Cycle 18). This would bring our 
total Lx/SFR inline with the prediction from Lehmer et al. (2019), 
within the errors. (For comparison, Lehmer et al. 2019 found a 
median intrinsic column density Мн = 2 x 10?! for their brighter 
sources, and applied that value to convert count rates to luminosities 
for the fainter sources). Even the choice of a higher but constant Мн 
would be an over-simplification (which is why we avoided it here): 
galaxies with a high sSFR contain more cold gas and should have an 
absorption column density above the median value (thus requiring a 
larger luminosity correction in Figure 11), while galaxies with low 
sSFR may be better approximated with Ny close to the line-of-sight 
Galactic value. 

Moreover, the relation between Lx, SFR and M. may be 
different for the galaxies in the Virgo sample and those in the 
Lehmer et al. (2019) sample. There are only 5 Virgo spirals 
(NGC 4254, NGC 4321, NGC 4450, NGC 4536, NGC 4569) in the 
sample of 38 galaxies used by Lehmer et al. (2019) to fit their scal- 
ing parameters. Even for those 5 galaxies, we cannot directly com- 
pare our and their adopted values of Lx, SFR and M., because 
Lehmer et al. (2019) use only the inner part of the D25 of those 
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Figure 12. Unabsorbed 0.3-10 keV luminosity of the most luminous ULXs 
with individual spectral modelling (Table 3) as a function of their spectral 
hardness (power-law photon index). Blue datapoints correspond to ULXs in 
the hard ultraluminous state (Г < 2); red datapoints correspond to the soft 
ultraluminous state (Г > 2). We omitted from the plot the unusual thermal 
source (young SN candidate) for which we cannot define a photon index 
(Section 4.2.2). 


spirals (about half of the projected 225 area). More specifically, 
Lehmer et al. (2019) consider only the portion of a galaxy inside the 
Ks = 20 mag (агсѕес) 2 surface brightness isophote (Jarrett et al. 
2003), while we have used the traditional D25 definition of galactic 
sizes. This may introduce a small bias, for example if the outskirts of 
spiral disks have a lower number of X-ray binaries and lower Lx for 
a given stellar mass or SFR, compared with the denser inner regions. 
Furthermore, we used different proxies for SFR and different prox- 
ies for stellar mass, compared with Lehmer et al. (2019). Finally, we 
used a different method for determining the total point-source lumi- 
nosity of a galaxy. We estimated Lx from a spectral analysis (Section 
3.1). Instead, Lehmer et al. (2019) estimate Lx by integrating the 
X-ray binary luminosity functions. Considering all these method- 
ological differences, the small discrepancy between our empirical 
estimates and the predicted scaling relations of Lehmer et al. (2019) 
is not surprising. A more detailed discussion of this issue is beyond 
the scope of this work. 


42 ULX population 
4.2.1 Spectral modelling for the ULXs 


In the first step of our ULX search (Section 3.2), we had assumed 
a fixed spectral model (power-law with photon index Г = 1.8) and 
Galactic line-of-sight column density to each host galaxy (typically, 
Мн = 2-3 x 1020 ст”; values from HIAPI Collaboration et al. 
2016 and Kalberla et al. 2005). Under those two assumptions, we 
identified 85 off-nuclear sources that exceed a 0.3-10 keV unab- 
sorbed luminosity of 1039 erg s^! in at least one observation. We 
report their positions, observed 0.5-8 keV fluxes and unabsorbed 
0.3-10 keV luminosities in Table 2, including the fluxes and lumi- 
nosities in the observations in which the same sources were below 
the ULX threshold. Table 2 provides a useful working list of inter- 
esting point-like sources; however, all luminosity values are likely 
to be under-estimated, because we have not yet accounted for the 
loss of X-ray photons due to intrinsic absorption in the host galax- 
ies and around the compact objects. We will now assess the effect 
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of intrinsic absorption, and then improve on our ULX catalogue by 
using individual spectral modelling at least for the brightest sources. 

We obtained reliable fits with free photon index and column 
density for all the sources with unabsorbed luminosities 24 x 1039 
erg 87! (Table 3). Instead, the majority of sources with Lx ~ 1-4 
х1039 erg s7! do not have enough counts below 1 keV to constrain 
Ny and Г (and therefore the unabsorbed luminosity) independently. 
This is particularly the case for sources observed in later Chandra 
cycles, with a much degraded ACIS sensitivity in the soft band. 
Thus, we limit our spectral analysis to the most luminous ULXs. 
At the signal-to-noise available for the majority of our sources, 
an absorbed power-law model provides an adequate fit, with plau- 
sible photon indices 1 < Г x 3. The only exception is a ULX 
(CXOU J122541.70+071338.9) in the disk plane of the perfectly 
edge-on galaxy IC 3322A, which appears dominated by thermal 
plasma emission, as we discuss later (Section 4.2.2). 

Based on our spectral modelling results, we found at least 25 
ULXs that reach a 0.3-10 keV luminosity of z4 x 10°9 erg s^! at 
least in one observation (Table 3); 9 of them reached a luminosity 
of 1049 erg 571, In addition, a 10th (І.Х exceeded 10*° erg sl, 
in NGC 4480, if its Tully-Fisher distance of 42 Mpc is correct; 
however, in that case, NGC 4480 would most likely not be a member 
of the Virgo cluster. Thus, we omitted that source from subsequent 
plots and modelling of the cumulative and differential luminosity 
distributions. For consistency, we also subtracted the SFR and stellar 
mass of NGC 4480 when we calculated the predicted contribution 
of HMXBs and LMXBs, and subtracted the projected area of its 
D325 when we calculated the predicted contribution of foreground 
and background sources. Moreover, Table 3 includes three sources 
whose identification as stellar-mass accreting compact objects is 
less clear: the peculiar thermal-plasma source in IC 3322A (Section 
4.2.2), and two sources that may be the nuclei of smaller galaxies 
behind or near NGC 4492 and NGC 4567 (Section 4.2.3). 

We examined the hardness distribution as a function of un- 
absorbed luminosity (Figure 12), for 28 individual observations (as 
listed in Table 3) of 25 different ULXs above a luminosity of 4 x 1039 
erg 571. We omitted from that plot the four above-mentioned sources 
with less clear identification, in the О25 of IC 3322A, NGC 4480, 
NGC 4492 and NGC 4567. In 15 out of 28 cases, the sources can 
be classified in the hard ultraluminous state (using the classification 
scheme of Sutton et al. 2013), with a photon index Г « 2 (which 
means that the spectral luminosity ELg rises with energy, in the 
Chandra band), and the other 13 are in the soft ultraluminous state, 
with Г > 2. The hardness-luminosity distribution (Figure 12) is also 
consistent with the alternative ULX classification of Pintore et al. 
(2014), into group 1 (harder) and group 2 (softer) sources. In both 
ULX classification schemes, the difference between between harder 
and softer sources is mostly attributed to the scattering optical depth 
of the super-critical disk outflow along our line of sight. In addition, 
a harder ULX spectrum may be the signature of a neutron star rather 
than BH accretor (Pintore et al. 2017; Walton et al. 2018). 

In our sample, there is no statistical difference (at least below 
1040 erg 871) between the luminosity of hard and soft sources, 
although the four most luminous sources (with Lx » 1.5 x 1040 
erg 571) are all in the hard state. The median photon index is ~1.8, 
which also justifies our initial assumption of the photon index. The 
median total absorption column density is 3.0 x10?! cm^?, that is 
a factor of 10 higher than what we had assumed in our preliminary 
luminosity estimates based on srcflux results (Col. 6 in Table 2). 

To correct (at least on a statistical level) the likely under- 
estimation of the preliminary luminosity estimates for the sources 
that are too faint for individual spectral modelling, we re-estimated 
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Table 2. Host-galaxy X-ray binary luminosities and ULX census 


Galaxy LXxRB.0.5-8 R.A.(J2000) Dec.(J2000) ОХ Бо 5-8 ULX І 03-10 ULX Lo 3-10 ObsID Obs. Date Detector Exp. Time Notes 
(1039 еге sl) (10-14 erg cm? 5 (1039 erg sl) (1039 erg 871) (ks) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (1) (12) 
МСС4178 5.4[4.6-6.1]* 12 12 44.51 +10 51 13.6 14.713.7-15.8 3.95[3.67—4.21 12748 2011-02-19 ACIS-S 36.3 
NGC 4192 6.4[4.5-8.8] 12 13 48.92 +14 55 01.1 11.2[9.6-12.8] 3.21[2.74—5.57 19390 2017-11-04 ACIS-S 14.9 
NGC 4197 10.3[8.0-12.7] 12 14 40.29 +05 49 00.0 2.62[1.69-3.82 2.76(1.78-4.03 19420 2017-07-26 ACIS-S 8.0 
12 14 41.31 +05 49 03.9 2.99[2.00—4.25 3.15[2.10—4.48 19420 2017-07-26 ACIS-S 8.0 
NGC 4206 1.59[0.53-2.75] 12 15 17.48 +13 00 15.5 2.721.79-3.93 1.46(0.96-2.11 19431 2017-02-06 ACIS-S 8.0 
NGC 4212 7.3[6.2-9:3] 12 15 35.53 +13 53 49.1 2.76[2.06-3.60 1.19[0.92-1.61 19395 2017-02-14 ACIS-S 14.9 
12 15 39.59 +13 53 42.2 13.0[11.2-14.7 5.77[5.00-6.58 а 19395 2017-02-14 ACIS-S 14.9 
NGC 4216 14.0[11.5-17.0] 12 15 50.05 +13 07 13.0 6.69[4.60—9.34 2.47[1.70-3.44 3.2[2.2-4.5 19391 2017-07-24 ACIS-S 9.5 
12 15 50.12 +13 06 18.3 7.07[5.46-8.96 2.61[2.01-3.31 4.8[2.8-21.0] 19391 2017-07-24 ACIS-S 9.5 
NGC 4222 2.7[1.0—5.3] 12 16 20.55 +13 18 06.5 2.70[1.73-3.95 1.88[1.21-2.76 24[1.6-3.5 19432 2018-05-10 ACIS-S 8.0 
NGC 4237 1.4[0.7-8.1] 12 17 09.86 +15 18 53.8 1.88[1.14-2.91 1.20[0.73-1.86 1.6[1.0-2.5 19404 2017-02-12 ACIS-S 8.0 
NGC 4254 16.6[15.5-17.7] 12 18 51.87 +14 24 21.6 3.23[2.17-4.59 1.03[0.69-1.46 1.7[1.1-2.4 7863 2007-11-21 ACIS-S 5.1 
3.17[2.65-3.71 1.01[0.84—1.18 1.3[1.1-1.5 17462 2015-02-16 ACIS-I 44.5 
1218 56.10 +14 24 19.3 46.0[41.6—50.5 14.6[13.3-16.1 21.6[19.0-24.8] 7863 2007-11-21 ACIS-S 54 
16.2(15.0-17.3 7.15[6.64—7.71 6.6[5.9—7.4 7462 2015-02-16 ACIS-I 44.5 
NGC 4276 «L5 - - - - 19427 2017-07-30 ACIS-S 52 
NGC 4293 4.9[3.4-6.4 = = = = 9409 2018-04-03 ACIS-S 8.0 
NGC 4298 4.0[2.0-6.6 12 21 36.10 +14 35 45.1 2.37[1.68-3.23 0.91[0.65,1.24] 1.2[0.9-1.7 9392 2018-04-09 ACIS-S 42 
3.51[2.38-4.96 1.35[0.91—1.90 9397 2018-04-09 ACIS-S 7.8 
NGC 4299 2.2[1.4-3.3 - - - - 7834 2007-11-18 ACIS-S 54 
NGC 4302 1.74[1.06—2.66] 1221 42.04 +14 3654.5 2.46[1.53-3.68 1.10[0.69—1.66 1.2[0.7-2.1 19392 2018-04-09 ACIS-S 42 
3.00[1.59—4.97 1.34[0.72-2.24 9397 2018-04-09 ACIS-S 7.8 
NGC 4303 6.9[6.0—7.9 12 22 01.24 +04 29 37.9 4.96[4.36-5.56 1.13[0.99-1.26 2149 2001-08-07 ACIS-S 28.0 
4.08[3.30-4.86 0.93[0.75-1.11 7550 2014-11-16 ACIS-S 9.8 
NGC 4307 ZAT - - - - 19405 2018-05-14 ACIS-S 47 
NGC 4312 0.32[0.12-0.55] - - - - 7083 2005-11-07 ACIS-S 1.9 
9414 2017-02-15 ACIS-S 8.0 
NGC 4313 0.9[0.1—1.9 - - - - - 9416 2018-04-14 ACIS-S 8.0 
NGC 4316 3.9[2.8—5.4 12 22 39.54 +09 20 21.0 3.44[2.61—4.43 3.60[2.73—4.63 5.4[3.0-22.0] 9400 2018-04-20 ACIS-S 4.9 
NGC 4321 12.0[10.9-13.1] 12 22 46.20 +15 48 50.0 1.45[1.16-1.74 0.57(0.46-0.69 0.9[0.7-1.1] 6727 2006-02-18 ACIS-S 37.9 
0.97[0.62—1.43 0.39[0.24—0.57 0.6[0.4—0.9] 9121 2008-04-20 ACIS-S 4.9 
2.73[2.11-3.47 1.08[0.83-1.36 1.7[1.3-2.1] 2696 2011-02-24 ACIS-S 4.9 
1.271.07-1.48 0.50[0.42-0.58 0.7[0.6—0.8] 4230 2012-02-16 ACIS-S 7.9 
1.35[0.72-2.25 0.53[0.28-0.88 0.7[0.4—1.1] 23140 2020-02-15 ACIS-S 0.0 
1.87[1.11—2.88 0.74[0.44—1.13 0.9[0.5—1.4] 23141 2020-03-13 ACIS-S 0.0 
12 22 47.37 +15 49 11.5 <0.07 <0.03 <0.05 6727 2006-02-18 ACIS-S 37.9 
5.93[4.94—6.93 2.33[1.95-2.73 3.6[3.0-4.2] 9121 2008-04-20 ACIS-S 4.9 
«0.23 «0.09 «0.1 12696 2011-02-24 ACIS-S 4.9 
«0.08 «0.03 «0.05 14230 2012-02-16 ACIS-S 7.9 
«0.3 «0.1 «0.1 23140 2020-02-15 ACIS-S 0.0 
«0.3 «0.1 «0.1 23141 2020-03-13 ACIS-S 0.0 
12 22 51.29 +15 51 23.3 <0.14 <0.06 <0.1 6727 2006-02-18 ACIS-S 37.9 
<0.03 <0.02 <0.03 9121 2008-04-20 ACIS-S 49 
0.55|0.25-0.99 0.22[0.10—0.39 0.3[0.1—0.5 12696 2011-02-24 ACIS-S 4.9 
4.75[4.31-5.19 1.87[1.70-2.04 2.6[2.4-2.9 14230 2012-02-16 ACIS-S 7.9 
1222 53.71 +15 48 46.2 «0.11 «0.05 «0.08 6727 2006-02-18 ACIS-S 37.9 
3.53[2.82-4.25 1.39[1.11—1.67 2.1[1.7-2.6 9121 2008-04-20 ACIS-S 4.9 
«0.19 «0.07 «0.1 12696 2011-02-24 ACIS-S 4.9 
«0.05 «0.02 «0.03 14230 2012-02-16 ACIS-S 7.9 
0.98[0.45-1.78 0.39[0.18-0.72 0.55[0.28-0.95] 23140 2020-02-15 ACIS-S 0.0 
0.34[0.05-0.88 0.13[0.02-0.34 0.12[0.08-0.35] 23141 2020-03-13 ACIS-S 0.0 
12 22 54.14 +15 49 12.6 4.02[3.55—4.49 1.58[1.40-1.77 2.4[2.1-2.7 6727 2006-02-18 ACIS-S 37.9 
3.8[12.4-15.3 5.30[5.88-5.99 14.4[8.1-31.1] 9121 2008-04-20 ACIS-S 49 
4.16[3.36-4.96 1.63[1.32-1.96 2.5[2.1-3.0 12696 2011-02-24 ACIS-S 4.9 
3.00[2.69—3.31 1.18[1.06—1.30 1.8(1.6-2.0 14230 2012-02-16 ACIS-S 7.9 
6.88[5.26-8.50 2.70[2.07-3.35 3.5[2.7-4.4 23140 2020-02-15 ACIS-S 0.0 
4.1[11.8-16.3 5.53[4.65-6.41 14.0[6.5—44.6] 23141 2020-03-13 ACIS-S 0.0 
12 22 54.26 +15 49 44.4 .87[1.37-2.49 0.74(0.54-0.98 1.0[0.8-1.4 6727 2006-02-18 ACIS-S 379 
2.12[1.57-2.78 0.83(0.62-1.09 1.3[1.0-1.7 9121 2008-04-20 ACIS-S 4.9 
2.43[1.55-3.56 0.96(0.61-1.40 1.3[0.8-1.9 12696 2011-02-24 ACIS-S 4.9 
2.47[2.09-2.85 0.97(0.82-1.12 1.3[1.1-1.5 14230 2012-02-16 ACIS-S 7.9 
2.57|1.69-3.71 1.01[0.67-1.46 1.3[0.9-1.9 23140 2020-02-15 ACIS-S 0.0 
.69[0.99—2.64 0.66[0.39-1.04 0.9[0.5—1.4 23141 2020-03-13 ACIS-S 0.0 
12 22 54.78 +15 49 16.3 6.79[6.15—7.43 2.67[2.42-2.92 4.2[3.8-4.6 6727 2006-02-18 ACIS-S 37.9 
1.1[9.72-12.4 4.34[3.82-4.88 6.9[5.4—10.4] 9121 2008-04-20 ACIS-S 4.9 
4.07[3.22-4.92 1.61[1.26-1.94 12696 2011-02-24 ACIS-S 49 
5.85[5.42-6.28 2.30[2.13-2.48 14230 2012-02-16 ACIS-S 7.9 
6.33[4.86-8.06 2.49[1.91-3.18 23140 2020-02-15 ACIS-S 0.0 
4.55[3.33-6.02 1.75[1.31-2.37 23141 2020-03-13 ACIS-S 0.0 
12 22 54.83 +15 49 18.4 3.04[2.60-3.49 1.20[1.02-1.37 6727 2006-02-18 ACIS-S 37.9 
4.53[3.63—5.43 1.79[1.39-2.08 9121 2008-04-20 ACIS-S 4.9 
5.11[4.19-6.03 2.01[1.65-2.38 12696 2011-02-24 ACIS-S 4.9 
3.63[3.28-3.98 1.43[1.29-1.57 14230 2012-02-16 ACIS-S 7.9 
4.14[2.98—5.54 1.62[1.18-2.18 23140 2020-02-15 ACIS-S 0.0 
2.90[1.94—4.13 1.14[0.76-1.62 23141 2020-03-13 ACIS-S 0.0 
12 22 55.69 +15 49 25.2 1.84[1.49-2.19 0.72[0.58-0.86 6727 2006-02-18 ACIS-S 37.9 
2.58[1.97-3.20 1.0110.77-1.26 9121 2008-04-20 ACIS-S 14.9 
2.88[2.21-3.57 1.13[0.87-1.40 12696 2011-02-24 ACIS-S 14.9 
1.93[1.62-2.25 0.76[0.64—0.88 14230 2012-02-16 ACIS-S 7.9 
2.49[1.63-3.59 0.98|0.65-1.42 23140 2020-02-15 ACIS-S 10.0 
1.91[1.16-2.92 0.75[0.46—1.15 23141 2020-03-13 ACIS-S 10.0 
12 22 57.47 +15 49 16.5 <0.07 <0.03 6727 2006-02-18 ACIS-S 37.9 
<0.34 <0.13 9121 2008-04-20 ACIS-S 14.9 
<0.13 <0.05 12696 2011-02-24 ACIS-S 14.9 
3.93[3.59—4.28 1.55[1.42-1.69 14230 2012-02-16 ACIS-S 7.9 
«0.75 «0.25 23140 2020-02-15 ACIS-S 10.0 
«0.46 «0.20 23141 2020-03-13 ACIS-S 10.0 
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Table 2 — continued Host-galaxy X-ray binary luminosities and ULX census 


Galaxy LXRB,0.5-8 
(1039 erg 571 ) 
(1) (2) 
NGC 4321 12.0[10.9-13.1] 
NGC 4330 1.62[0.46—3.66] 
NGC 4343 0.5[0.3-1.6] 
NGC 4356 «1.0 
NGC 4380 2.7[2.0-4.0] 
ІС 3322A 6.4[4.6—8.5] 
NGC 4388 4.6[3.7-5.6] 
NGC 4390 0.3[0.1-1.1 
IC 3322 0.6[0.16—1.6] 
NGC 4394 1.0[0.5-2.0 
NGC 4396 3.8[2.9-4.7 
NGC 4402 2.3[1.7-3.4 
NGC 4405 1.0[0.6—2.2 
NGC4411A 1.49[0.88-2.25] 
NGC 4412 3.6[1.8—6.6 
NGC 4411B 2.4[1.7-3.6 
МСС 4407 2.7[0.9—4.9 
МСС 4416 <1 
МСС 4419 4.3[3.2-6.5 
NGC 4424 0.9[0.5-1.5 
NGC 4429 4.1[2.0-6.3 
NGC 4430 2.6[1.9-3.7 
NGC 4438 3.9[2.9-5.1 
NGC 4445 1.0[0.4—2.2] 
NGC 4450 8.4[5.6-11.3] 
NGC 4451 1.3[0.7-2.1] 
1C 3392 «0.6 
NGC 4457 0.7[0.4—0.9] 
NGC 4469 1.1[0.6-1.6] 
NGC 4470 1.0[0.8-1.2] 
NGC 4480 11.0[6.8-16.6] 
NGC 4492 3.7[2.6—5.2] 


NGC 4496A 17.1[12.0-23.3] 


NGC 4498 0.9 
NGC 4501 12.7[ 
МСС 4519 1.9 
NGC 4522 


0.5-14] 
9.5-15.9]* 


0.9-2.9] 
«0.9 
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R.A.(J2000) 


(3) 


12 22 57.52 


12 22 58.64 


12 22 58.68 


12 23 13.26 


12 25 22.21 
12 25 41.70 
12 25 42.25 
12 25 38.98 


12 25 39.14 


12 25 54.06 


12 26 01.54 
12 26 01.81 


12 26 27.86 
12 26 31.62 


12 26 35.16 
12 26 48.76 


12 26 49.48 


12 27 27.09 


12 30 27.86 
12 3028.14 
12 30 57.83 


12 30 58.88 
12 31 38.08 


12 32 00.39 
12 32 00.94 
12 32 06.19 


Dec.(J2000) 


(4) 


+15 48 47.5 


+15 51 24.0 


+15 47 51.9 


ES 


12126.0 


001 03.4 


+07 13 38.9 
+07 13 19.8 


РЕ 


$ 


2 40 00.8 


2 40 09.4 


+12 40 04.2 


+15 39 41.4 
+15 39 41.8 


+08 51 45.03 
+08 52 28.1 


+03 58 15.0 
+08 53 38.1 


+08 52 18.9 


+06 15 53.9 


+09 26 08.4 
+17 04 26.8 


+07 49 32.3 


+04 14 59.9 
+04 14 40.9 
+08 04 35.2 


+08 04 52.2 


+03 56 42.2 


+14 24 42.4 
+1425 02.6 
+1423 21.5 


ОХ Fọ.5-8 


(10:14 еге cm? sl) 


«0.073 


(5) 


«0.08 
«0.13 
«0.13 
«0.06 
3.55-6.27 
3.74—6.54 
«0.13 
«0.02 
«0.1 
2.92-3.65 
«0.38 
2.79-3.60 
2.98-4.42 
1.84-3.06 
2.86-3.47 
1.21-3.01 
0.93-2.58 
2.51-5.05 


2.16-3.83 
4.26-6.67 
071-1.84 
2.37-3.41 


3.44-4.67 
1.13-1.88 
1.29-2.07 
2.48-3.48 


2.97-5.66 
10.1-14.6 


2.29-4.67 
2.55-5.50 
0.51-1.90 
2.49-5.02 
0.36-1.62 
0.96-2.70 
1.98-3.91 
1.85-4.10 
2.79-5.03 


5.53[4.04-7.35] 


2.86 
3.68) 


3.11 
4.23 
2.57| 
1.86 
2.10) 
2.27 
2.78) 
1.14 
4.02 
0.83 
7.01 
4.18 
2.54| 
35:5 


8.05 
9.16 
2.36 


2.12-3.76 
2.58-5.04 


2.46-3.78 
3.12-5.35 
1.84-3.44 
1.35-2.47 
1.57-2.71 
1.51-3.24 
2.11-3.56 
5.89-1.96 
2.89-5.42 
0.33-1.65 
6.05-7.99 
2.98-5.67. 
1.94-3.15 
29.3-41.8 


7.09-9.02 
8.15-10.2 
1.76-3.05 


ULX Lo.3-10 


(1039 erg sl) 


2.10[1.53-2.79] 


48 
30 


(6) 


<0.03 
<0.05 
<0.05 
<0.03 
1.40-2.47 
1.47-2.56 
<0.05 
<0.01 
<0.05 
1.15-1.44 
«0.15 
1.13-1.46 
1.18-1.74 
0.73-1.20 
1.12-1.36 
0.48-1.18 
0.37-1.01 
1.42-2.86 


121-2.16 
4.01-6.39 
0.68-1.78 
1.19-1.71 
<0.037) 

1.72-2.35 
0.57-0.94 
0.66-1.06 
1.28-1.79 


0.77-1.46 
2.58-3.76 


0.98-2.00 
1.09-2.34 
0.21-0.81 
1.06-2.13 
0.70-3.20; 
0.72-2.05 
1.50-2.01 
1.40-3.10 
2.08-3.74 


1.09-1.94 
0.92-1.79 


Chandra survey of Virgo spirals 


ОХ LC 


(103 


23 
25 
0.7 
19 
24 
2.0 
28 
34 

55] 


2.9 


1.9 
1.8 


0.3-10 


9 erg sl) 
(7) 


<0.05 
<0.08 
<0.08 
<0.05 
1.8-3.2 
1.9-3.3 
<0.08 
<0.02 
<0.08 
1.6-2.0 
<0.20 
1.8-2.4 
2.0-2.8 
1.2-1.9 


1.6-3.2 
1.7-3.6 
0.3-1.3 
1.3-2.6 
11-45 
11-32 
2.0-3.9 
22-48 
2.8-29.3] 


[2.13.8] 


ObsID 


(8) 


6727 
9121 
2696 
4230 
23140 
23141 
6727 
9121 
2696 
4230 
23141 
6727 
9121 
2696 
4230 
23140 
23141 
9419 
4687 
7129 
2955 
19435 
9406 
9401 
9401 
1619 
2291 
1619 
12291 
1619 
12291 
19425 
19424 
7864 
19417 
19417 
15149 
19410 
7837 
7840 
7837 
19429 
19418 
7840 
19429 
7840 
19429 
19434 
19436 
19394 
19408 
19430 
19426 
2883 
8042 
21376 
23189 
23037 
23200 
19433 
19399 
19412 
19421 
3150 
19415 
12978 
15756 
15760 
16260 
16261 
16262 
21647 
19423 
19423 
7845 
15759 
7845 
15759 
16995 
19422 
2922 
2922 
2922 
19411 
19428 


Obs. Date 


(9) 


2006-02-18 
2008-04-20 
2011-02-24 
2012-02-16 
2020-02-15 
2020-03-13 
2006-02-18 
2008-04-20 
2011-02-24 
2012-02-16 
2020-03-13 
2006-02-18 
2008-04-20 
2011-02-24 
2012-02-16 
2020-02-15 
2020-03-13 
2018-04-16 
2005-02-11 
2006-07-16 
2011-02-17 
2017-04-16 
2018-05-11 
2018-04-24 
2018-04-24 
2001-06-08 
2011-12-07 
2001-06-08 
2011-12-07 
2001-06-08 
2011-12-07 
2018-04-14 
2018-04-22 
2007-11-11 
2018-05-25 
2018-05-25 
2013-07-24 
2018-04-09 
2007-11-13 
2007-02-13 
2007-11-13 
2017-02-09 
2017-11-23 
2007-02-13 
2017-02-09 
2007-02-13 
2017-02-09 
2018-05-10 
2018-05-06 
2016-12-24 
2017-04-17 
2018-05-10 
2018-05-05 
2002-01-29 
2008-02-11 
2020-03-20 
2020-03-20 
2020-03-27 
2020-03-28 
2018-05-13 
2017-04-06 
2018-05-13 
2018-04-11 


2007-02-22 
2014-04-25 
2007-02-22 


2002-12-09 
2002-12-09 
2002-12-09 
2018-05-05 
2018-05-05 


Detector 


(10) 


ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-I 
ACIS-I 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-I 
ACIS-S 
ACIS-I 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 


Exp. Time 


(ks) 
(11) 


37.9 
14.9 
14.9 
79 
10.0 
10.0 
37.9 
14.9 
14.9 
79 
10.0 
37.9 
14.9 
14.9 
7:9. 
10.0 
10.0 
8.0 
38.3 
4.6 
72.8 
14.9 
147 
13.9 
13.9 
20.0 
27.6 
20.0 
27.6 
20.0 
27.6 
15.4 
14.9 
5.1 
8.1 
8.1 
148.1 
8.0 
5.1 
4.9 
зі 
9.9 
8.0 
4.9 
9.9 
4.9 
9.9 
7.8 
7.8 
9.9 
14.9 
7.8 
7.8 
25.1 
4.9 
29.7 
19.8 
19.8 
25.7 
14.9 
8.0 
14.7 
7.6 
38.9 
14.9 
19.8 
32.1 
29.4 
24.7 
22.8 
24.7 
29.7 
8.0 
8.0 
4.9 
29.7 
4.9 
29.7 
2.7 
8.1 
17.9 
17.9 
17.9 
8.5 
78 


15 


Notes 


(12) 


>>>>>> 


mmo 
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Table 2 — continued Host-galaxy X-ray binary luminosities and ULX census 


Galaxy LXRB,0.5-8 R.A.(J2000) ^ Dee.(J2000) ULX Fo 5-8 
(1039 erg 871) (10:14 еге cm 2 sl) 
(1) (2) (3) (4) (5) 
NGC 4527 7.5[5.6—9.6] 2 34 10.94 +02 39 25.2 24.7[21.6-27.9 
0.43[0.13-0.95 
2 34 11.43 +02 39 28.9 18.0[15.4—20.7 
3.11[2.20-4.23 
1234 11.76 +02 39 57.7 «0.34 
3.94[2.91-5.18 
2 34 14.87 +02 40 05.3 <0.34 
5.28[4.08—6.72 
NGC 4532 11.9[10.6-13.4] 2 34 18.89 +06 28 15.4 20.0[17.1—23.0 
2 34 20.34 +06 27 37.5 11.4[9.22-13.7 
NGC 4535 3.1[1.2-7.5] 2 34 29.94 +08 10 35.4 7.68[6.35—9.02 
NGC 4536 12.0[10.8-13.1] 2 34 26.16 +02 11 22.3 9.31[7.89-10.7 
2 34 26.83 +02 11 18.8 3.08[2.30-4.02 
2 34 31.81 +02 11 21.9 3.99[3.05—4.93 
12 34 32.56 +02 10 54.5 7.89[6.57-9.21 
12 34 36.93 +02 08 51.4 3.1[11.3-15.0 
12 34 38.81 +02 10 28.5 5.57[4.41-6.74 
NGC 4548 22[1.7-3.4] 2 35 31.33 +14 30 48.8 4.72[3.37—6.37 
NGC 4567 9.6[7.2-13.4] 12 36 33.75 +11 16 10.6 4.96[3.62—6.59 
NGC 4568 12.1[9.0-15.1] 12 36 32.32 +11 13 06.3 3.18[2.10-4.54 
3.411.74-5.82 
12 36 33.04 +11 1405.6 2.1[9.72-14.4 
7.3[13.6-20.8 
12 36 33.22 +11 14 37.1 9.87[7.82-11.9 
.90[0.92—3.42 
12 36 34.48 +11 1440.2 2.20[1.54—3.03 
.86[1.08-2.93 
NGC 4569 7.4[6.3-8.5 1236 53.70 +13 11 54.0 4.93[3.03-7.50 
2.53[2.16—2.89 
МСС 4571 2.0[0.3-5.6 - - - 
NGC 4579 28.3[26.5-30.2]* 12 37 40.31 +11 47 27.8 28.9[27.8-30.1 
12 37 43.21 +11 4901.5 3.34[2.53-4.15 
12 37 53.88 +11 50 20.4 9.47[8.14—10.8 
МСС 4580 1.210.6-2.1 - - - 
NGC 4606 1.2[0.4-2.1 - - - 
NGC 4607 3.2[2.4-4.7 1241 12.61 +11 53 50.2 1.49[0.75-2.60 
7.98[6.12-9.85 
NGC 4639 5.2[3.8-7.8 12 42 51.20 +13 14 40.3 3.52[1.82-6.05 
4.70[3.66—5.75 
NGC 4647 6.9[5.1—8.8 - - - 
NGC 4654 3.9[3.0-4.9] - - - 
NGC 4689 5.7[4.0-7.4] 12 47 47.71 +13 46 17.3 7.93[6.09-9.79] 
NGC 4698 9.0[6.8—11.5] 12 48 25.87 +08 30 20.8 14.7[13.8-15.7] 
NGC 4713 2.0[1.5-3.4] - - - 


ULX Lo.3-10 ULX Lo 3-10 ObsID Obs. Date Detector Exp. Time Notes 
(1039 erg 871) (1039 erg 871) (ks) 
(6) (7) (8) (9) (10) (11) (12) 
7.4(6.5-8.4| 21.4[18.1-25.7] 4017 2003-03-09 ACIS-S 4.9 
0.13[0.04—0.29 0.18[0.06—0.40] 9386 2016-12-18 ACIS-S 9.9 
5.4[4.6-6.2] 13.3[11.1-16.1] 4017 2003-03-09 ACIS-S 4.9 
0.94[0.66—1.27 1.3[0.9-1.8 19386 2016-12-18 ACIS-S 9.9 
«0.10 «0.17 4017 2003-03-09 ACIS-S 4.9 
1.18[0.87-1.56 1.6[1.2-2.1 9386 2016-12-18 ACIS-S 9.9 
<0.10 <0.17 4017 2003-03-09 ACIS-S 4.9 
1.58[1.22-2.01 22[1.7-2.7 9386 2016-12-18 ACIS-S 9.9 
5.13[4.38-5.91 7.3(5.9-10.5| 9407 2018-05-15 ACIS-S 8.0 
2.93[2.36-3.51 7.3[3.3-27.9] 9407 2018-05-15 ACIS-S 8.0 
(2.98[2.47-3.50]) (4.0[3.3-4.7]) 9388 2017-04-2 ACIS-S 14.9 F 
3.15[2.67-3.64 5.1[4.1-7.3 9387 2017-07-01 ACIS-S 14.9 
1.08[0.80—1.40 1.4[1.0-1.9 9387 2017-07-0 ACIS-S 14.9 
1.35[1.04-1.66 1.8[1.4-2.3 9387 2017-07-0 ACIS-S 14.9 
2.67[2.23-3.12 3.7[3.1—4.3 19387 2017-07-0 ACIS-S 14.9 
4.44[3.80—5.06 9.8[5.2-24.0] 9387 2017-07-0 ACIS-S 14.9 
1.89[1.49-2.28 2.5[2.0-3.0 19387 2017-07-0 ACIS-S 14.9 
1.77[1.25-2.38 2.4[1.7-3.2 19402 2018-04-1 ACIS-S 7.6 
3.87[2.83-5.14 (4.9[3.4—7.4]) 19396 2018-05-10 ACIS-S 7.8 E 
1.4410.95-2.05 1.6[1.1-2.1 19389 2018-03-0 ACIS-S 13.7 
1.58[0.81—2.72 2.0[1.0-3.5 19396 2018-05-10 ACIS-S 7.8 
5.43[4.38-6.50 6.7[5.2-8.6 19389 2018-03-0 ACIS-S 13:7 
7.81[6.14—9.43 9.0[7.0—13.6] 19396 2018-05-10 ACIS-S 7.8 
4.45[3.43—5.36 4.9[3.9-6.2 19389 2018-03-0 ACIS-S 13.7 
0.86(0.42-1.55 1.6[1.0-2.3 19396 2018-05-10 ACIS-S 7.8 
1.03[0.72-1.14 1.4[0.9-1.9 19389 2018-03-0 ACIS-S 13:7 
0.87[0.51—1.37 1.2[0.7-1.8 19396 2018-05-10 ACIS-S 7.8 
1.00[0.61—1.52 1.8[1.1-2.7 405 2000-02-17 ACIS-S 17 
0.51[0.44—0.59 0.82[0.70-0.94] 5911 2005-11-13 ACIS-S 39.1 
- - 7858 2008-02-14 ACIS-S 47 
12.6[12.1-13.1 17.1[16.2-18.1] 807 2000-05-02 ACIS-S 33.9 
1.41[1.07-1.82 2.4[2.1-2.7 807 2000-05-02 ACIS-S 33.9 
4.12[3.53-4.74 4.0[2.8—4.7 807 2000-05-02 ACIS-S 33.9 
- - 10413 2017-12-0 ACIS-S 8.0 
- - 19437 2017-04-22 ACIS-S 77 
0.89[0.44—1.54 1.2[0.6-2.2 19437 2017-04-22 ACIS-S 77 
4.74[3.64—5.87 10.5[5.5-50.6] 19403 2018-05-09 ACIS-S 8.0 
2.58(1.33-4.45 4.8[2.5-16.8] 408 2000-02-05 ACIS-S 14 
3.45[2.68—4.22 4.5[3.5-5.5 19398 2018-04-14 ACIS-S 14.7 
- - 785 2000-04-20 ACIS-S 38.1 
8182 2007-01-30 ACIS-S 52.3 
8507 2007-02-0 ACIS-S 17.5 
12975 2011-08-08 ACIS-S 84.9 
12976 2011-02-24 ACIS-S 101.0 
14328 2011-08-12 ACIS-S 14.0 
- - 3325 2002-05-03 ACIS-S 4.9 
19393 2016-12-14 ACIS-S 9.9 
2.99[2.29-3.68] 4.7[3.6-5.8] 7865 2007-05-07 ACIS-S 4.9 
NA NA 3008 2002-06-16 ACIS-S 29.7 G 
- - 4019 2003-01-28 ACIS-S 4.9 


Col. (2): total unabsorbed 0.5-8 keV luminosity of the point-source (X-ray binaries) population (resolved plus unresolved) in each galaxy, estimated with a 
power-law model fit. The range of values in brackets represents the 90% confidence range. When a galaxy was observed multiple times, the luminosity was 
determined from a stack of all the observations (combined with the epicspeccombine task). An asterisk (on 3 galaxies) denotes that they are larger than the 
detector field of view in one direction. See Section 3.1 for details. Cols. (3),(4): coordinates of each ULX, defined as each source that exceeded a 0.3-10 keV 
luminosity of 10?? erg s7! in at least one of the observations. Col. (5): absorbed 0.5-8 keV flux (with 90% confidence range) of each ULX in each 
observation, measured with srcflux; here, we assumed a power-law model, with fixed photon index Г = 1.8 and only Galactic line-of-sight column density. 
Col. (6): unabsorbed 0.3-10 keV luminosity (with 90% confidence range) of each ULX in each observation (including those in which they were in a faint 
state), assuming Г = 1.8 and Galactic line-of-sight column density. The conversion from fluxes to luminosities is based on the median redshift-independent 
NED distances. The error range in the luminosity includes only the flux measurement uncertainty, not the systematic uncertainty in the distance to each 
galaxy. Col. (7): “corrected” 0.3-10 keV ULX luminosity (with 90% confidence range), based on individual spectral models for the most luminous sources, 
and, for all the others, adopting the median values of photon index and column density inferred from the modelled subsample (Г = 1.8 and total 


Ny = 3.0 x 10?! , respectively). We regard these luminosities as our best estimate, and use them for further analysis. Col. (12): Notes. A = SN 1979C; В = 
likely background galaxy at z « 0.156; C = likely foreground star; D = NED Tully-Fisher distance possibly overestimated; Е = candidate nuclear source of a 


small (satellite?) galaxy; F = likely quasar; G = known BL Lac at z « 0.43. 


all their fluxes again with srcflux, and Г = 1.8, but this time with 
the assumption of Ng iy = 3.0 х 10?! cm, for all the sources 
not already included in Table 3. In summary, the list of "corrected" 
0.3-10 keV luminosities provided in Table 2, Col. (7), is our best es- 
timate, combining the spectral fitting results for the most luminous 
sources and the revised srcflux estimates for all the others. On av- 
erage, the corrected unabsorbed luminosities for sources observed 
in earlier Chandra cycles (before the soft-energy degradation) are 
a factor of 1.6-1.7 times higher than estimated when only Galac- 
tic line-of-sight absorption is assumed. In Cycle 18, the difference 


between Мн = 3.0 x 102! cm"? and Ny = 3.0 x 1029 ст”? was 
reduced to a factor of «1.3. 


4.2.2 Candidate Type Пп SN in IC 3322A 


CXOU J122541.70+071338.9, detected in IC 3322A (distance of 
x25 Mpc) on 2018 April 24, is the most luminous ULX in the whole 
sample. Its X-ray spectrum (Table 3) suggests multi-temperature 
thermal plasma emission and an intrinsic column density of «2 
х1022 erg s~!, for an intrinsic luminosity of «6 x10% erg 571. A 
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spectrum dominated by X-ray thermal plasma at such a high lumi- 
nosity is very unusual for a ULX. The X-ray spectrum is consistent 
with that of a young SN; more specifically, of a Type IIn, which 
explodes in a denser circumstellar medium. Based on the Super- 
nova X-ray Database online catalog? (Ross & Dwarkadas 2017), 
there are at least five core-collpase SNe (SN 1998Z, SN 1995N, SN 
2005ip, SN 2005kd, SN 2006jd and SN 2010]]), detected at X-ray 
luminosities of a few 1040 erg 571, even up to a few 10*! erg s^! , 
for several years after the explosion. See also Chandra et al. (2012, 
2015) and Chandra (2018) for detailed discussions of the X-ray and 
multi-band properties of Type IIn SNe. The possibility that some of 
the brightest ULXs may be young SNe in a dense environment was 
discussed in other cases, for example for the brightest ULX in the 
Cartwheel galaxy (Pizzolato et al. 2010). 

No core-collapse SNe have ever been reported in IC 3322A, 
but that is not strong evidence against our suggested interpreta- 
tion. The peak absolute optical brightness of Type IIn SNe is 
broadly distributed between about —17 mag and -21 mag, with 
a median value around -19 mag (Nyholm et al. 2020). The dis- 
tance modulus of IC 3322A is «32 mag. To that, we may add 
10 mag of visual extinction, based on the best-fitting value of 
Мн = 2x102 cm"? (Predehl & Schmitt 1995; Güver & Özel 2009; 
Willingale et al. 2013). Thus, it is likely that the presumed SN was 
fainter than 20 mag in the optical, and would have been easily missed 
by large-area sky searches. 

We used the Open SN Catalog? catalogs (Guillochon et al. 
2017) to check whether any of the point-like sources (even those 
fainter than 10°? erg 571, not listed in Table 2) may correspond to 
a historical SN. There are 39 optically-identified SNe reported to 
occur in the 75 sample galaxies prior to at least one of the Chandra 
observations used for this study. Only three of them are detected 
in X-rays (Table 4). The most luminous of them is SN 1979C (a 
Type II L) in NGC 4321, for which we infer an X-ray luminosity of 
жб x10°8 erg s^! when its spectrum (stacked over the six Chandra 
observations of this galaxy) is fitted with a two-temperature thermal 


plasma model. The best-fitting temperatures are КТ = Dag 


keV and KT» = (2775 keV, and the best-fitting intrinsic column 


density is (5857) x 10?! cm?, 


4.2.3 Likely interlopers from other projected galaxies 


Two bright sources (Table 3) are associated with the nuclei of (ap- 
parently) small galaxies projected within the D»s of a larger spiral 
(NGC 4492 and NGC 4567). We found no measurements for the 
redshifts of the two small companion galaxies. If the host galaxies 
are dwarf satellites at approximately the same distance as NGC 4492 
and NGC 4567, respectively, the two X-ray sources would have a 
luminosity of several times 1039 еге 87! (Table 3). However, they 
would not be classified as ULXs in the most common definition 
of this class, because they are nuclear sources. They would still 
be interesting as candidate intermediate-mass BHs in the nuclei 
of dwarf satellites. Instead, the bright source projected inside the 
D»s of NGC 4698 (Tables 2, 3), known in the literature as XMMU 
J124825.9+083020 or WISE J124825.85+083020.4, is a BL Lac 
at redshift 0.43, confirmed by optical spectroscopy (Foschini et al. 
2002). 
Another 


curious source (Table 3) is CXOU 


? https://kronos.uchicago.edu/snax/ 
10 һирѕ://ѕпе.ѕрасе/ 
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Figure 13. Cumulative LF іп the 0.3-10 keV band. See Col. (7), Table 2 for 
the data values. The datapoints include a subsample of ULXs modelled with 
free Ny and Г (Table 3), while for the fainter sources we fixed Г- 1.8 and 
Nu,tot = 3-0 X 10?! cm~?. The red histogram represents the luminosities 
inferred from the earliest exposure available for each galaxy with multiple 
observations; the green histogram shows the luminosities inferred from 
the latest available exposures for each galaxy; the blue histogram is for 
the longest available exposures. The small excess at the higher end of the 
earliest-look LF is the result of small-number statistics (namely, two ULXs 
with Lx > 10*° erg 87! in NGC 4527, detected only in the first observation 
of that galaxy). For graphical purposes only, the most luminous source 
(candidate SN) in the Virgo sample (Section 4.2.2) is falling outside the 
plotted region, at Lx ж 6 x 10* erg 571. The black dashed line is the 
predicted contribution from the cosmic background plus foreground stars 
(Lehmer et al. 2012). The black dotted line shows that the lower end of the 
cumulative LF is consistent with a power-law of index —1, with a downturn 
at the high-luminosity end. 
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Figure 14. Differential LF in the 0.3-10 keV band. For each of the four 
luminosity bins, of equal width Alog (Lx /erg s~!) = 0.3, the colours are 


defined as in Figure 13 (red, green and blue for the earliest, latest and 
longest observation, respectively). The relative position of the three coloured 
datapoints in each bin is offset by a very small amount along the X axis 
for visual clarity. The predicted background and foreground contributions 
have already been subtracted from each luminosity bin. The dashed black 
line is an indicative power-law model, with the form dN /d(log L) œ 15 
with y = —0.80; values of y = -0.8 + 0.2 are consistent with the observed 
datapoints. For each luminosity bin, a black line marks the value predicted for 
that bin by the LF models of Lehmer et al. (2019) (LMXBs plus HMXBs). 
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Figure 15. Comparison between the empirical cumulative LF from our 
survey, and the model components from Lehmer et al. (2019). Here, we 
have plotted all functions in the 0.5-8 keV band, to facilitate the comparison. 
Conversions between 0.3-10 keV and 0.5-8 keV bands were done according 
to the best-fitting spectral models of the sources. Red histogram: cumulative 
ULX LF for the Virgo spiral sample (here we chose the distribution derived 
from the earliest observations). Dashed blue line: predicted contribution 
from HMXBs, for a star formation rate of 70M yr! (similar to the SFR 
estimated by B15). Dashed red line: predicted contribution from LMXBs in 
spiral disks, for a stellar mass of 1.5 x 1012 Mo. Dotted black line: predicted 
contribution from background galaxies, AGN, and foreground stars. Solid 
green line: sum of the previous three terms (HMXBs, LMXBs and bg/fg 
sources), which represents the expected observed luminosity function from 
our Virgo sample. Dotted magenta line: predicted contribution from LMXBs 
also for a stellar mass of 1.5 х 102 Мо, but calibrated on the younger 
population of elliptical galaxies (Zhang et al. 2012). For graphical purposes 
only, the most luminous source (candidate SN) in the Virgo sample (Section 
4.2.2) is falling outside the plotted region, at Lx = б x 109 erg 571. 


J122538.98+124000.8, projected inside the D25 of NGC 4388, 
but also projected onto the spiral arm (not the nucleus) of an 
anonymous background galaxy at redshift z ~ 0.156, seen behind 
the Virgo galaxy. If this transient source belongs to NGC 4388, it 
is a run-of-the-mill ULX with Lx « 1.5 x 1039 erg s instead, 
if it belongs to the background spiral, it is a hyperluminous 
X-ray source with Lx ж 2.6 x 1042 erg s^! (in the 2001 June 
8 observation), a strong intermediate-mass BH candidate. In a 
subsequent observation of the same galaxy a decade later (2011 
December 7), the source was no longer detected (Table 2), which 
implies a decline of a factor of at least 40. 


4.3 Properties of the ULX luminosity function 


LFs are ideally defined and quantified as instantaneous snapshots of 
an X-ray population. Collecting and including all sources (including 
transient ones) detected in several successive observations of the 
same galaxy leads to an over-estimate of the luminosity distribution. 
For this reason, when plotting the 0.3- 10 keV cumulative LF (Figure 
13), we included only sources from one observation per galaxy. 
We compared the distribution obtained by selecting the earliest 
observation for each galaxy, the latest observation, and the longest 
observation. The three alternative histograms are consistent with 
each other at the low luminosity end. At the high end, the small 
excess of the earliest-look LF is the result of small-number statistics: 
two ULXs with Lx > 109 erg 67! in NGC 4527, detected only in 
the first observation of that galaxy. 


We used the log N-log S curves in the Chandra Deep Field 
study of Lehmer et al. (2012) to estimate the number of contami- 
nating sources (background galaxies, AGN, quasars, and foreground 
stars) in our survey. To do this, we calculated the expected number 
of such sources projected inside the D25 of each of our sample 
galaxies, above a flux corresponding to a 0.3-10 keV luminosity of 
1039 erg s^! for that galaxy (at its median NED distance, given in 
Table 1, Column 8). For most of the galaxies, the expected number is 
<1. Adding those numbers for all the galaxies, we estimate a total of 
~15 contaminating sources above 1039 еге s7! (Figure 13). In fact, 
the total number of 15 is likely to be an over-estimate, because it as- 
sumes that the count rates and fluxes of background AGN suffer no 
attenuation from photoelectric absorption passing through the halo 
or disk of the Virgo galaxies. With a similar method, we estimate 
the expected contamination at other luminosity levels; for example, 
only about 2 contaminating sources are expected above ~5 x10°9 
erg 871. (A few candidate background/foreground sources have al- 
ready been flagged in Table 2 and Table 3). 

We divided the energy range between 1.4 x 1039 erg 871 and 
2.2 x 10% erg 87! (0.3-10 keV band) into four equally spaced bins 
(in log scale), and computed the contribution of the three alterna- 
tive LFs (earliest, latest and longest) in each bin. We then subtracted 
the expected background and foreground contributions in each bin 
from each of three observed distributions. This gives three alter- 
native realizations of the differential LF of the point-like sources, 
corresponding to three different snapshots in time (Figure 14). As 
expected, there is a lot of scatter at the high-luminosity end, due to 
small-number statistics. A power-law with an index y = -0.80--0.20 
(defined as dN /d (log L) сс Д7) is a good approximation of the slope 
of the differential LF. This is intermediate (as expected, given the 
mix of the two populations) between the slopes of ~ —0.6 for a 
pure HMXB population, and ~ —1.3 for LMXBs in that energy 
range (Lehmer et al. 2019). The normalization is also consistent 
(Figure 14) with the spiral galaxy models of Lehmer et al. (2019), 
for a total SFR «70 Мо yr! (B15 value) and a total stellar mass 
М, tot = 1.5 x 10? Мо. This does not imply that the SFR of 100 
Мо уг”! inferred from the WISE photometry is wrong. It simply 
means that the SFR proxies (far-UV GALEX plus 24um Spitzer 
luminosity) used by Lehmer et al. (2019) to calibrate the normal- 
ization of the HMXB population provide SFR values more similar 
to the B15 values than to the WISE-derived values. 

Our empirical differential LF falls slightly below the model 
predictions at the low-luminosity end (Figure 14). This suggests 
that we are missing a few sources. This is expected. Sources with 
a 0.3-10 keV luminosity of ~1.5 x10°9 erg s^! and an absorbing 
column Ng z a few times 1021 cm"? have a substantial chance 
of missing our first selection cut (Table 2, Col. 6), depending on 
their photon index and on the Chandra observing cycle. For this 
work, we are not particularly concerned by this issue: ULXs in the 
luminosity range ~(1-1.5) x 1039 erg s^! are the upper end of the 
ordinary X-ray binary population; we are more interested in the 
high-luminosity end of the ULX population. 

The total 0.3-10 keV luminosity of the resolved D25 popula- 
tion above 1.4x 10?? erg 5—1 is -(3.5-3.7) x10! erg s^! (depending 
on the choice of observation for galaxies with multiple exposures). 
About 10% of this luminosity comes from background AGN pro- 
jected inside the D25 regions. An extrapolation of the differential 
LF (Figure 14) with a power-law slope of y = —0.8 + 0.2 down to 
1036 erg s~! predicts an additional contribution of «2-3 x10*! erg 
871 from unresolved point sources. This value is not very sensitive 
to the choice of the lower limit, as most of the integrated emission 
comes from the high-luminosity sources. Thus, our best estimate for 
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Table 3. Spectral parameters of the most luminous ULXs in Virgo (Lx > 4 x 10°? erg 871) 


Host Galaxy 


NGC 4178 


NGC 4197 


NGC 4212 


NGC 4216 


NGC 4254 


NGC 4316 


NGC 4321 


ІС 3322A 


NGC 4396 


NGC 4411B 


NGC 4496А 


NGC 4501 


NGC 4527 


NGC 4532 


NGC 4536 


NGC 4568 


NGC 4579 


NGC 4607 


NGC 4480 


NGC 4492 


NGC 4567 


NGC 4388 


NGC 4698 


8 


В.А (12000) 


12 12 44.51 
2 14 40.29 
2 1441.31 
2 15 39.59 
2 15 50.12 


12 18 56.10 


2 22 39.54 


2 22 54.14 


2 22 54.78 
225 41.70 
22601.81 
2 26 49.48 
231 38.08 
2 32 00.39 
2 32 00.94 
2 34 10.94 
2 3411.43 
2 34 18.89 
2 34 20.34 
2 34 26.16 
2 34 36.93 


2 36 33.04 


2 36 33.22 
2 37 40.31 


2 37 53.88 


241 12.61 


12 30 28.14 


12 25 38.98 


12 48 25.87 


Dec.(12000) 


+10 51 13.6 
+05 49 00.0 
+05 49 03.9 
+13 53 42.2 
+13 06 18.3 


+14 24 19.3 


+09 20 21.0 


+15 49 12.6 


+15 49 16.3 
+07 13 38.9 
+15 39 41.8 
+08 52 18.9 
+03 56 42.2 


+14 24424 


+14 25 02.6 


+02 39 25.2 


+02 39 28.9 


+06 28 15.4 


+06 27 37.5 


+02 11 22.3 


+02 08 51.4 


+11 1405.6 


+11 1437.1 


+11 47 27.8 


+11 50 20.4 


+11 53 50.2 


+04 14 40.9 


+08 04 35.2 


+11 16 10.6 


+12 40 00.8 


+08 30 20.8 


ObsID 


12748 
9420 
9420 
9395 
19391 
7863 


7462 


9400 


19429 
16995 
2922 
2922 
4017 
4017 
9407 
9407 
9387 
9387 
9389 
9396 
9389 
807 
807 


9403 


19423 


15759 


19396 


1619 


3008 


Obs. Date 


2011-02-19 


2017-07-26 


2017-07-26 


2017-02-14 


2017-07-24 


2007-11-21 


2015-02-16 


2018-04-20 


2008-04-20 


2020-03-13 


2008-04-20 


2018-04-24 


2018-05-25 


2017-02-09 


2015-07-26 


2002-12-09 


2002-12-09 


2003-03-09 


2003-03-09 


2018-05-15 


2018-05-15 


2017-07-01 


2017-07-01 


2018-03-01 


2018-05-10 


2016-12-06 


2014-04-25 


2018-05-10 


2001-06-08 


2002-06-16 
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Мн сәй Мин г ма C-stat/dof Fesa 
(1022 cm 2) (1022 cm?) (10714 erg cm? 571) 
Most likely ULXs 
+0.09 0.21 +0.9 42.2 
0.026 0.24" o 09 1.477051 3.3703 206.9/235 19.7250 
+1.21 +2.5 +5.8 +5.8 
0.015 0.237753 1.324 0.6754 11.5/14 4376 
+0.8 +0.44 +2.6 
0.015 <0.8 0.8706 0.307916 30.8/17 4.6750 
+0.25 +0.46 +1.8 +3.1 
0.026 0.197019 1.59 0743 2.7710 92.3/104 14.075: 
+1.13 +1.15 +10.2 +3.0 
0.028 0.817577 2.0175 96 3.07; 26.8/36 TA Y 
+0.10 +0.28 +3.1 +10.3 
0.026 0. 107 0:09 1.357037 8.8» 137.4/170 64.398 
+0.16 +0.25 +1.3 +1.9 
0.026 0.2714 1.627 0774 3.679 185.4/237 17.3717 
Д +0.53 +1.13 32.7 +1.3 
0.014 0.2724 2.44* 04 13% 22.9/34 2.9109 
+0.12 +0.61 +3.6 +1.2 
0.018 0.2875 14 3.467035 7.3753 101.5/112 8.811 
+0.33 +0.81 +9.1 +3.0 
0.018 0.321027 3.1 10:72 8471 54.7/75 10.97; 
40.16 40.48 42.4 42.3 
0.018 0.3170 14 2.28 0 14 4.5715 87.2/115 10.6773 
0.017 20051-21 (NA-thermal source)? 30.9/46 319* US 
+0.61 +0.92 +17.3 +2.8 
0.020 0.73048 2.84 0:80 9.67: 55.0/61 10.2775 
+0.72 +1.26 jr5.2 *1.7 
0.013 0.547 0:51 2.399 08 1.9713 31.6/28 3.5115 
+0.52 43.5 +11.8 
0.023 < 0.20 1.207030 4.8714 61.4/80 44.27 193 
+0.14 +0.77 42.7 +0.7 
0.027 0.311012 3.8370 66 4.7716 91.5/91 41707 
+0.12 +0.42 +1.5 +1.8 
0.027 0.269 14 2.20 039 3.3710 130.2/127 8.9715 
40.46 40.50 7.6 +22.1 
0.020 0.84 0735 1.237046 7.9 6 99.1/126 55.3 19.9 
40.35 40.49 4.2 +8.8 
0.020 0.557 0797 1.227 0:45 4.9771 71.0/99 36.4 75 
+0.37 +0.53 +4.8 +4.9 
0.016 0.407930 1.737049 5.4754 89.3/99 21.7741 
+0.43 +0.93 +10.1 +3.5 
0.016 0.437 0733 2.97 0.79 7.0 38 37.7/55 8.7170 
+0.49 +0.62 +2.8 +3.0 
0.016 0.527 0:39 1.57 0.51 2.475 86.9/89 1.4753 
+0.28 +0.64 +5.9 +2.3 
0.016 0.2975: 3.0075 7. 10 64.9/92 10.2724 
+0.37 +0.54 +0.87 +3.2 
0.024 0.057 0705 0.827534 0.887530 68.6/75 13.0755 
+0.72 +0.74 44.2 44.6 
0.024 0.47 0 47 1.3679 63 2.5715 52.0/66 15.5756 
+0.21 +0.36 42.4 
0.024 «0.16 1.207051 1.097059 71.4/78 10.2775 
+0.03 +0.12 +0.9 42.5 
0.031 0.1170 03 1.727015 74^ 8 275.4/300 35-7 34 
+0.11 +0.51 +1.1 +2.0 
0.031 0.1070 09 2.017073 2.195 92.5/92 7.8776 
+0.88 +1.12 +12.7 +3.0 
0.025 0.8275) 62 2.3795 4.6734 36.9/40 7.631 
Likely ULX but very uncertain distance 
+0.58 +1.2 +2.32 32.3 
0.021 0.177017 1.8778 0.9675 55 26.3/19 4.7435 
Nuclei of apparently small galaxies: either dwarf satellites in Virgo, or unrelated background galaxies at unknown redshift 
+0.54 +0.57 +3.4 +1.8 
0.014 0.9672: 1.877533 3.115 89.6/108 8.715 
+0.94 +4.2 +1.57 
0.024 < 0.67 12274 2.5 14 29.6/28 0.707 39 
Either a ULX in Virgo, or an off-nuclear hyperluminous X-ray source іп a background galaxy at z = 0.155 
40.28 +0.11 +1.0 
0.026 <0.08 1.07057 0.57 0:09 51.8/71 З.А в 
Background BL Lac at z ғ 0.43 
40.04 40.21 40.7 41.5 
[0.016] 0.07* 0 04 2.05*0 1o 3.8705 164.5/210 14.2773 


Р. observed fluxes in the 0.5-8 keV band. 


шо 


©: if located at the distance of NGC 4492. 


“қ 


8: if located inside NGC 4388. 


h 
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: if located at the distance of NGC 4567. 


: units of 1075 photons keV~! cm"? s^! at 1 keV. 


0.32 


г isotropic de-absorbed luminosities in the 0.3-10 keV band, defined as 4rd? times the de-absorbed 0.3-10 keV fluxes. 


: thermal-plasma spectrum, fitted with 2-temperature apec model; kT, = 0.641047 keV; kT = 2.0 keV (frozen). 


: if projected onto a spiral arm of a background galaxy (z ж 0.155: Alam et al. 2015) seen through the D25 of NGC 4388. 


С: 
103-10 


(1039 erg 571 


+0.6 
6.0752 


45.9 
5.2724 


3.8 
4.9717 


+1.8 
TA 


+16.2 
4.8 50 


43.2 
21.676 

+0.8 
6.6757 


+16.6 
5.41161 


+16.7 
14.41.5 


+30.6 
14.0304 


43.5 
6.97; 


+69.8 
62.2288 


331.0 
94%; 


+23.8 
5.5755 


+5.2 
16.67€ 


+22.5 
12.9422. 


+2.1 
6.0 1 


44.3 
21.473 
adi 
7.3+3-2 


-1.4 


420.6 
73740 


12.2 
5.1750 


+14.2 
9.8746 


+1.9 
6.7215 


+4.6 
9.0750 


+1.3 
4.910 


+1.0 
Tarne 


+1.2 
4.0703 


+40.1 
10.5005 


+22.4 
13.0722; 


40.518 
(1.67031) 


h 
4743 
(2652143) 


(1.4104 x 10°) 
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Table 4. Historical SNe observed in the sample 


Host Galaxy 


NGC 4178 
NGC 4254 


NGC 4302 
NGC 4303 


NGC 4316 
NGC 4321 


NGC 4402 
NGC 4411b 
NGC 4419 


NGC 4424 


NGC 4451 


NGC 4496a 


NGC 4501 
NGC 4527 


NGC 4536 


NGC 4568 


NGC 4579 


NGC 4639 
NGC 4647 


SN 


19631 
1967H 
1972Q 
19861 


2014L 
1986E 
1926A 
19611 
1964F 
1999gn 
2006оу 
2008in 
2014dt 
2003bk 
901B 
1914A 
959E 
979C 
2006X 


202001 


976B 
1992ad 
1984A 
2012cc 
1895A* 
2012cg 
1985G 
1960F 
1999с1 
1915A 
1991T 
2004gn 
1981B 
1990B 
2004сс 
1988А 
1989М 
1990N 
1979А 


R.A.(J2000) 


12 12 45.61 
12 18 55.12 
12 18 50.64 
12 18 52.04 


12 18 48.68 
1221 41.21 
12 21 54.15 
12 22 00.45 
12 21 53.02 
12 21 57.02 
12 21 55.30 


1221 57:57. 
12 22 42.14 
12 22 47.60 
12 22 57.00 
12 22 58.91 
12 22 58.63 
12 22 53.99 
12 22 54.96 
12 26 11.91 
12 26 49.59 
12 26 55.71 
12 26 56.81 
12 27 16.90 
12 27 12.83 
12 28 40.40 
12 31 42.06 
12 31 56.01 
12 34 10.95 
12 34 10.17 
12 34 12.10 
12 34 29.58 
12 36 33.83 
12 36 34.40 
12 37 43.54 
12 37 40.72 
12 42 56.70 
12 43 29.14 


Dec.(J2000) 


+10 5131.0 
+14 24 40.3 
+14 26 36.3 
+14 24 44.1 


+14 24 43.5 


+14 37 55.0 


+04 28 18.5 
+09 20 01.2 
+15 49 25.0 
+15 47 30.0 
+15 49 01.3 
+15 4751.7 
+15 48 33.1 
+15 49 25.0 
+13 06 44.8 
+08 52 38.2 
+15 03 17.7 
+15 02 45.5 
+09 25 05.0 
+09 25 13.2 
+09 15 40.0 


+03 56 47.9 
+14 25 35.3 
+02 39 03.5 
+02 39 56.4 
+02 39 34.4 
+02 11 59.3 
+11 1429.8 
+11 14 32.8 
+11 48 19.4 
+11 49 26.0 
+13 15 23.7 
+11 35 27.1 


Туре 


Peak Date 


963-05-05 
967-07-01 
972-12-17 
986-05-17 
2014-01-26 
986-04-13 
926-05-09 
961-06-03 
964-06-13 


999-12-17 
2006-11-24 
2008-12-26 
2014-10-29 
2003-02-28 
901-02-10 
914-03-02 
959-02-21 
979-04-15 
2006-02-19 
2020-01-07 
976-04-05 
992-07-01 
984-01-17 
2012-04-29 
895-03-16 
2012-05-17 
985-03-17 
960-04-20 


999-06-15 
915-03-20 
991-04-28 
2004-12-01 
98 1-03-09 


990-01-18 
2004-06-10 
1988-01-19 
1989-06-28 
1990-07-10 
1979-01-15 


a 
Po s-8 


(10716 erg cm? 871) 
<7 
<8 
<8 


<8 


(20:2 ) 13 


+26 
3875 


< 10 
< 10 
< 10 
< 10 
x55, 
< 15 
< 15 
<37 
<10 
<7 

<7 
16020 
<7 
« 1004 
«17 
< 30 


< 30 


+50, 
50730 


<25 
«25 
«20 
« 80 
«15 
< 20 
< 20 
< 30 
< 30 
< 30 
< 30 
<5 
<5 
< 20 


<3 


b 
Lo 5-8 


(1037 erg 871) 
<1.5 
<2.0 
<2.0 
<2.0 


<10 


4; observed fluxes in the 0.5-8 keV band. To estimate the non-detection upper limits, we assumed а thermal-plasma apec model, with fixed temperature 
kT = 1 keV and line-of-sight Galactic Ny. For SN 2014L, SN 1986E, SN 2012cc, apec temperatures and column densities were left as free parameters in the 


fit, with Cash statistics. For SN 1979C, a two-temperature apec model was used (Section 4.2.2). 


Р. de-absorbed luminosities in the 0.5-8 keV band, based on the models mentioned above. 
©: likely due to extended thermal emission 


4; located in a region of strong diffuse emission 


©: first historical SN discovered outside the Local Group 
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the background-subtracted point-source 0.3-10 keV luminosity in 
the Virgo sample is -5-6 х104! erg s^, corresponding to ~4.5-5 
x104! erg s7! in the 0.5-8 keV band. 

First, we compare this extrapolated value with the integrated 
background-subtracted point-source galaxy luminosities estimated 
in Section 4.1. There, we had obtained a luminosity of «3.6 x 104 
erg s^! in the 0.5-8 keV band. However, for many of the galaxies 
with moderate or low emission (or short exposures) we had simply 
assumed a Galactic line-of-sight absorption in the conversion from 
count rates to fluxes and luminosities. Consequently, that value must 
be considered a lower limit. We have already seen (Section 4.2.1) 
that increasing the column density from 3 x 1029 cm? to 3 x 10?! 
cm ? increases the inferred luminosity by 20 to 50 per cent (for a 
given count rate). Therefore, we conclude that the total luminosity 
estimates from the integration of the D25 source regions and from 
the (extrapolated) point-source LF are consistent with each other. 
We had already noted (Section 4.1) that the value of ~3.6 x104! 
erg s7! appeared under-estimated for the given mass and SFR of 
the sample, and suggested that this under-estimation particularly 
affected galaxies with dusty disks, seen at high inclination. 

Second, we compare the measured point-source luminosity 
with the model predictions (Lehmer et al. 2019), in the energy range 
where we have better constraints (Lx = 1.2 x 1039 еге 571 in the 
0.5-8 keV band). Direct sum of the ULX luminosities (minus the 
predicted AGN and foreground star fraction) gives a net 0.5-8 keV 
luminosity of «2.7 x10*! erg 571, but largely dependent on the 
small-number statistics of the most luminous sources. Over the 
same energy range, for an SFR of 70 Мо уг”! and a stellar mass 
of 1.5 x 1012 Мо, model predictions are 43.2 х104 erg 871, of 
which «2.5 x10! erg 87! from HMXBs and «0.7 x 104! erg s^! 
from LMXBs. (The number of HMXBs per unit luminosity interval 
starts to dominate over the number of LMXBs at Lx ~ 1.2 x 1039 
erg s^). As we noted earlier, the observed ULX luminosities are 
more consistent with the lower end of the published SFR values 
(B15 values) rather than the WISE-derived values. 

We analyzed the 0.5-8 keV cumulative luminosity distributions 
of observed sources and model predictions (Figure 15) to determine 
whether the small excess luminosity in the model of Lehmer et al. 
(2019) is significant. For sources below «4 x10°9 erg 571, the fidu- 
cial model follows our empirical distribution extremely well. Above 
that energy bin, the model is significantly above the observed LF. 
We suggest several alternative reasons for the discrepancy. The first 
possibility is that there is an imperfect match between the sector 
of the LF for which most of the sources have an individual spec- 
tral fit, and the sector where the luminosity was estimated with a 
fixed Ny and Г. Splicing the two sub-samples together may cre- 
ate a normalization jump between them. The second possibility is 
that our simple (power-law) spectral models are still underestimat- 
ing the luminosities of a few of the most luminous sources (more 
than for the fainter ones). For example, the most luminous ULXs 
might have a stronger thermal component below 1 keV, in addition 
to the power-law component, which was not included in the spectral 
models. More likely, our luminosity function may include a few 
ULXs with an intrinsic luminosity Lx > 4 х 1039 erg s^! but an 
absorption column > 1022 ст”? (е. g., іп dense star-forming regions 
or inside an edge-on galactic disk), which have been incorrectly 
placed in lower-luminosity bins. Finally, and perhaps more interest- 
ingly, the assumed shape of the model functions may not be correct. 
In Lehmer et al. (2019), both the LMXB and HMXB differential 
LFs in the ULX range were assumed to be unbroken power-laws 
(with different slopes) up to a cut-off at 5 x 1049 erg 871, While 
this may be justified for HMXBs, there is no sufficient evidence to 
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say that it is also a good model for LMXBs. In fact, in the ULX 
sample studied by Lehmer et al. (2019), there are no LMXBs above 
a luminosity of ~2 x10°9 erg s~!. If the LMXB distribution has 
a cut-off at Lx « a few x10?? erg 51, it would help explain why 
our observed ULX distribution is consistent with the LMXB plus 
HMXB model function at the low luminosity end but follows the 
HMXB-only function at the high luminosity end (Figure 15). 


4.4 ULXSsin early-type and late-type spiral galaxies 


Let us now compare the total number of ULXs found in our study 
with the numbers predicted by the models of Kovlakas et al. (2020) 
for spiral galaxies (valid from 50/а to Sm). The energy band used to 
define the ULX luminosity in their study is 0.5-8 keV. In the same 
band, we observed «60 ULXs, after removing the expected back- 
ground and foreground contributions. (We used the first observation 
of each galaxy for this estimate, but the number is approximately the 
same if we choose the longest observation or the latest observation.) 
Considering a mild effect of incompleteness around 1.0-1.3 x 1039 
erg 571, hinted at in Section 4.3, the corrected number of ULXs at 
any given time is probably ~65—70, but we keep the estimate of 60 
as a safe lower limit. The model of Kovlakas et al. (2020) predicts 


Мх © 0.45+0:06 х SFR 13.3725 X M, 15, where М, 12 is in units 


of 1012 Мо. For a meaningful comparison with our numbers, we 
need to consider the definition of SFR used by Kovlakas et al. (2020) 
for their calibration of the ULX scaling relation. In their case, the 
SFR was derived from the total infrared luminosity, primarily from 
the /RAS photometric measurements (calibration of Dale & Henou 
2002; Kennicutt & Evans 2012) supplemented for some galaxies by 
the 12-um апа 22-um WISE measurements. Therefore, for consis- 
tency, we use the WISE measurements of SFR for our galaxy sample 
rather than the B15 values (the latter being more biased in favour 
of UV- and Ha-based SFR measurements). Namely, we take a total 
SFR ~ 86-104 Mo уг”!, where the lower and higher values corre- 
sponds to the 22-um and 12-um estimates, respectively (Table 1). 
For this SFR range, the scaling relation of Kovlakas et al. (2020) 
predicts a best-fitting Ny_x = 44-52. If we consider also the er- 
ror range in the coefficients of the scaling relation, the predicted 
number of ULXs in our sample is ~36—58, lower than observed. 
Kovlakas et al. (2020) did point out that their ULX luminosities and 
hence ULX numbers are lower than those found in other surveys, 
including that of Swartz et al. (2011) (which is the most similar to 
the present work, for methodology). The difference was attributed 
to a different morphological composition of the samples, and a dif- 
ferent conversion procedure between Chandra/ACIS count rates, 
fluxes and luminosities (Appendix СІ in Kovlakas et al. 2020). 
More importantly, a comparison with the scaling of 
Kovlakas et al. (2020) is interesting because of their claimed strong 
difference between ULX rates in early-type spirals (SO/a to Sbc) and 


late-type spirals (Sc or later). For the former, Nur x = 0.1670 x 


SFR +11.2*32 x M, 12; for the latter, Мух = 0.98*0-1 х SFR. 
To test those relations, we split the Virgo spiral sample into two 
corresponding early-type and late-type subsamples, based on the 
NED classification (see also Table 1 and Figure 4). The early- 
type-spiral sample includes 45 galaxies, with a total stellar mass 
of «1.2 x10!2Mo and а WISE-derived SFR of 455-66 Мо yr-!; 
the late-type-spiral sample has 30 galaxies, with a stellar mass of 
«0.3 x10? Мо and SFR «30-38 Мо yr`!. For the parameters of 
the early-type-spiral sample, Kovlakas et al. (2020) predict a cen- 
tral value of 23-24 ULXs (516-31 when we include the uncer- 
tainty on the scaling coefficients); our observed number is 39. For 
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the late-type-spiral sample, the predicted central value is 429—37 
ULXs («24-42 with full uncertainty range); our observed number 
is 21. We conclude that the scaling models of Kovlakas et al. (2020) 
underpredict the number of ULXs in early-type spirals and over- 
predict those in late-type ones, at least when applied to the Virgo 
sample. 


5 DISCUSSION AND SUMMARY 


We are conducting the first large-scale Chandra/ACIS survey of the 
spiral galaxy population in the Virgo cluster. This project comple- 
ments and integrates the AMUSE-Virgo survey, which studied the 
early-type galaxy population of Virgo, a decade ago. We selected 
a sample of 75 galaxies, including all the largest and most active 
spirals («20 late-type galaxies with SFR >1 Mo), plus a selection 
of fainter ones. As detailed in Section 2.1, our sample is skewed to- 
wards brighter objects, in the sense that about 60 of our 75 galaxies 
have SFR z0.3 Mo yr-!. This was done to provide a contrast with 
the early-type galaxy sample (Gallo et al. 2008; Plotkin et al. 2014; 
Graham & Soria 2019) and to aid the discovery of luminous XRBs. 

The total observing time was about 1.95 Ms (combining new 
and archival observations), and for all galaxies in our sample we 
reached a point-source detection limit of «3 x10°8 erg 571 (or 
deeper, in a few cases). 

In this first paper, we have outlined the general properties of the 
sample: the distance, morphological type, stellar mass and SFR. We 
have then presented the X-ray population properties of the sample 
galaxies, in relation to their stellar mass (for which we estimate a 
total of «1.5 х1012 Мо) and SFR. For the SFR, we used the WISE 
12um апа 22um proxies and compared them with a collection of 
other SFR proxies (in particular, based on Ha and the far-UV) 
in the literature. Estimates vary between «70 Mo yr! and «100 
Mo yr !. Two-thirds of the galaxies can be described as active 
star-forming disks (dominated in X-rays by the HMXB population), 
while the others are early-type spirals, with low sSFR and an X-ray 
luminosity dominated by LMXBs. 

We identified a sample of 85 ULXs that exceeded a 0.5-8 
keV luminosity of 1039 erg 571 in at least one observation, and 
determined their coordinates. We fitted the spectral parameters for 
the most luminous sources (25 ULXs with Lx = 4х 1039 erg 
671) and used their median power-law index and column density to 
convert count rates to luminosities for the rest of the ULXs, when 
they were too faint for individual fitting. We showed that the ULX 
luminosity distribution from our Virgo study is broadly consistent 
with the LF models from Lehmer et al. (2019), for a SFR ~ 70Mo 
уг”! and a total stellar mass of «1.5 х1012 Мо. 

The model functions of Lehmer et al. (2019) predict a sig- 
nificant contribution (in number) from LMXBs: about half of 
the sources in the ~(1-1.5) x10?? erg s7! range, and about one 
third of our total number of ULXs. At first sight, this appears to 
contradict the estimate from Plotkin et al. (2014) (based on their 
AMUSE-Virgo study), of only about 1 old-population ULX per 
М. = 1.6 x 10!! Mo, all of them with Lx x 2 x 1039 erg s7!. (In 
our case, this would correspond to only about 10 old-population 
ULXs rather than about 30). In fact, there is no contradiction. The 
term “LMXB” is generally applied to a variety of X-ray binary popu- 
lations that are proportional to stellar masses, but the normalization, 
break and slope of their LFs are substantially different, depending 
on age and environment (Zhang et al. 2012). In the disk of spiral 
galaxies, especially for sSFR «ТІП уг, the number and lumi- 
nosity of bright LMXBs (above a few times 1037 erg 871) per unit 


mass is a few times higher than in old ellipticals. This is evident for 
example in the properties of the luminosity function in the disk of 
the nearby disk-dominated spiral M 83, where the LMXB compo- 
nent is at least as important as the HMXB component (Long et al. 
2014). Most of the stellar mass in the AMUSE-Virgo survey is in 
massive elliptical galaxies (Table 1 in Gallo et al. 2010). Therefore, 
the scaling relations of (Zhang et al. 2012) are suitable to model 
the LMXB population mapped by that survey. Instead, the LMXB 
population in our sample includes a mix of spheroidal (spiral bulge) 
population for the early-type spirals, and pure disk population. Thus, 
we expect it to follow the scaling relation derived by Lehmer et al. 
(2019), at least up to a few times 1029 еге sl, 


The contribution of LMXBs to the higher-luminosity ULXs 
remains an open question. From the models of Lehmer et al. (2019) 
we expect four old-population sources above Lx « 4x 1039 еге 57! 
(0.5-8 keV band), compared with about 16 HMXBs. Instead, our 
Virgo ULX LF is consistent with the HMXB component alone. A 
cut-off of the LMXB population at Lx = 4 х 1039 erg 571 is one of 
the possible explanations we have proposed for the slight mismatch. 
Some theoretical population synthesis models (Wiktorowicz et al. 
2017) do predict the formation of neutron star ULXs at least as 
luminous as 1040 erg s^! during phases of steadily super-critical 
mass transfer from low-mass donors (main sequence, Hertzsprung- 
gap stars, red giants, asymptotic giants, or white dwarfs), long after 
the end of star formation. The predicted number of old-population 
UL Xs does, however, drop rapidly with age, decreasing by a factor 
of 20 from 1 to 5 Gyr after a star formation burst, and by another 
order of magnitude from 5 to 10 Gyr (Wiktorowicz et al. 2017). 
In addition to those steady mass-transfer systems, old-population 
ULXs should include also transient systems with a large (outer 
radius z 1012 cm) accretion disk, subject to the thermal-viscous disk 
instability, which may exceed the Eddington limit during outburst 
(Hameury & Lasota 2020). 


In further work currently in preparation, we will illustrate and 
discuss the optical counterparts and stellar environments of the 
ULXs in this Virgo galaxy sample, based on the Next Generation 
Virgo Cluster Survey maps (Ferrarese et al. 2012) from the Canada- 
France-Hawaii Telescope. We will distinguish between younger and 
older populations, for the lower-luminosity and higher-luminosity 
subsets of ULXs, and for the early-type-spiral and late-type-spiral 
subsamples of galaxies. We will also discuss the metallicity effect 
on the high-luminosity end of the younger population. 


Individual spectral studies of the most interesting Chandra 
sources will also be presented in follow-up work. Here, we have 
simply mentioned a few highlights. The most luminous ULX in 
the sample (a source in IC3322A with Lx ~ 6x 10% erg s7!) 
may not be an accreting compact object, as its X-ray spectrum 
is more consistent with that of a young SN (which could have 
been missed by optical/IR searches). Apart from that mysterious 
source, the three most luminous ULXs (in NGC 4254, NGC 4496A 
and NGC 4579) all reached a luminosity Lx ~ 2 х 10% erg 571, 
which is the well-known characteristic threshold above which ULXs 
become much rarer (Grimm et al. 2003; Swartz et al. 2004, 2011; 
Mineo et al. 2012a). At least three Chandra sources seen inside 
the D25 of sample galaxies are likely to be the nuclear sources 
of smaller galaxies: further investigation is needed to determine 
whether those galaxies are also in Virgo (possibly satellites of their 
larger companions) or are a more distant background. 
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